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Preface 


This is a companion volume to Electromagnetism for Electronic Engineers (3™ edn.) (Ventus, 2009). 
It contains the worked examples, together with worked solutions to the end of chapter examples, 
which featured in the previous edition of the book. I have discovered and corrected a number of 
mistakes in the previous edition. 


I hope that students will find these 88 worked examples helpful in illustrating how the fundamental 
laws of electromagnetism can be applied to a range of problems. I have maintained the emphasis on 
examples which may be of practical value and on the assumptions and approximations which are 
needed. In many cases the purpose of the calculations is to find the circuit properties of a component 
so that the link between the complementary circuit and field descriptions of a problem are illustrated. 


Richard Carter 
Lancaster 2010 
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1. Electrostatics in free space 


1 Electrostatics in free space 


1.1 Introduction 


Electrostatic problems in free space involve finding the electric fields and the potential distributions 
of given arrangements of electrodes. Strictly speaking “free space’ means vacuum but the properties 
of air and other gases are usually indistinguishable from those of vacuum so it is permissible to 
include them in this section. The chief difference is that the breakdown voltage between electrodes 
depends upon the gas between them and upon its pressure. The calculation of capacitance between 


electrodes in free space is deferred until Chapter 2. 


The other problems included in this chapter involve the motion of charged particles (electrons and 
ions) in electric fields in vacuum. This topic remains important for certain specialised purposes 


including high power radio-frequency and microwave sources, particle accelerators, electron 


microscopes, mass spectrometers, ion implantation and electron beam welding and lithography. 


1.2 Summary of the methods available 


Note: This information is provided here for convenience. The equation numbers in the companion 
volume Electromagnetism for Electronic Engineers are indicated by square brackets. 


Symbol Signifies 
& (epsilon) The primary electric constant 


a) Electric charge 

q Electric line charge 

o (sigma) Surface charge density 

p (rho) Volume charge density 

E Electric field 

V Electric potential 

V (del) The vector differential operator 
T,X, Y,Z Unit vectors 


Inverse square law of force between charges in free space 


£00... 


F r 


Aner 


e Definition of the electric field of a charge in free space 


Units 

8.854 x 10°? F.m'! 
C 

Cai’ 

C.m” 

C.m° 

Vim" 

Vv 


[1.1] 


[1.2] 
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e Force acting on a charge placed in an electric field 
F=Q,E [1.3] 


e Gauss’ Theorem 
The flux of E out of any closed surface in free space is equal to the charge enclosed by the surface 
divided by €o. 


e The integral form of Gauss’ Theorem 
1 
E.dA =—||| pdv [1.3] 
| € N) 


e The differential form of Gauss’ Theorem 


i 6 6OCOE F 
ve-(° 2 he }- ie [1.9] 


ox Oy E 


0 

e Electrostatic potential difference 
B 

V,-V,=-| E-dl [1.13] 


e Calculation of electric field from the electrostatic potential 


E= die cane =-—grad V =-VV [1.22] 
Ox Oy Oz 


e Poisson’s equation 


2 2 2 
om 
ON a _ P 


vy / [1.24] 
Ox” Oy? Ozx” Es 
e Laplace’s equation 
gy OV OV OV 
VV= =0 [1.27] 


Sa 5 va 
Ox” Oy Oz 


e The Principle of Superposition and the method of images 
e The Principle of conservation of energy 
e The finite difference method 
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Example 1.1 


Find the force on an electron (charge -1.602 x 10°’ C) which is 1 nm from a perfectly conducting 
plane. What is the electric field acting on the electron? 


Solution 


Using the method of images the conducting plane is replaced by an image charge of +1.602 x 10'°C 
which is | nm behind the position of the conducting plane 


The force acting on the electron is found using the inverse square law [1.1] noting that the charges are 
2 nm apart. 


19,2 
) 


— OQ, _ (1.602 x10 
~ dee -12 _9\? 
0” 4nx8.854x10 x(2x10 


12 


F =-57.7x10 “~ N (1.1) 


Force is a vector quantity so a complete answer must specify its direction. The negative sign indicates 
that the electron is attracted to the image charge. The force is therefore acting towards the plane and at 
right angles to it. 


The electric field acting on the electron is found by substituting its charge and the force acting on it 
into [1.3] 


12 
19 


F -57.7x10— 
Q -1.602x10— 


1 


E= = 360 MV-m_ (1.2) 


The electric field is a vector quantity and the positive sign indicates that it is acting away from the 
plane. 


Example 1.2 


The surface charge density on a metal electrode is o. Use Gauss’ theorem to show that the electric 


field strength close to the surface is E = o/é, . 


Solution 


Consider a small element of area of the surface dA such that the surface around it can be considered to 
be a plane. The local charge density can be considered to be constant and, from symmetry 
considerations, the electric field must be normal to the conducting surface. Now construct a Gaussian 
surface dS, as shown in fig. 1.1, such that it encloses the element dA and has sides which are normal 
to the surface and top and bottom faces which are parallel to the surface. 
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Fig. 1.1 A Gaussian surface for calculating the electric field of a surface charge. 

Since E is parallel to the sides of dS the flux of E through the sides is zero. Also, because the electric 
field within a conducting material is zero when the charges are stationary, the flux of E through the 
bottom of dS is zero. The flux of E through the top of dS is 

d®=EdA (1.3) 
where E is the magnitude of E (since E is normal to the top of dS). The total charge enclosed by dS is 


dQ=adA (1.4) 


By Gauss’ theorem 


dD = ae (1.5) 
“0 
Substituting in (1.5) from (1.3) and (1.4) gives 
pa (1.6) 
E 


Note: Because a conducting surface is always an equipotential surface when the charges are stationary 
E must always be normal to it. If the surface is curved the electric field varies over it (1.6) shows that, 
locally, the charge density is always proportional to the electric field. 


Example 1.3 


Figure 1.2 right shows a charged wire which is equidistant from a pair of earthed conducting planes 
which are at right angles to each other. 

a) Where should image charges be placed in order to solve this problem by the method of images? 
b) What difference would it make if the planes were at 60° to each other? 

c) Could the method be used when the planes were at 50° to each other? 
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a 


d 


Fig. 1.2 A charged wire close to the intersection of two conducting planes 
Solution 
a) If Cartesian co-ordinates are used to describe the positions of the wire and of its images in the 


plane then the image line charges are —g at (- d, d) and (d, - d) and +¢ at (- d, - d) as shown in fig. 
1.3. 
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b) 


| 
eee eames 
| 


: 
o 


Fig. 1.3 Image charges for planes intersecting at 90° 


c) When the planes are at 60° to each other five image charges are equally spaced on a circle as 
shown in fig. 1.4. 


GF 
\ 7 ; 

! \ Toe \ 
alee ee pe 

\ aX 

cae a 
Y/i\ oO 
\. ie . 7 
SS, 


Fig. 1.4 Image charges for planes intersecting at 60° 


d) No. The method can only be used when the angle between the planes divides an even number of 
times into 360°. Thus it will work for planes at angles of 1/4, 1/6, 1/8, 1/10 of 360° and so on. 


Example 1.4 


A wire | mm in diameter is placed mid-way between two parallel conducting planes 10 mm apart. 
Given that the planes are earthed and the wire is at a potential of 100 V, find a set of image charges 
that will enable the electric field pattern to be calculated. 


Solution 


If we were to put just one image charge on either side of the wire the field pattern could be calculated 
by superimposing the fields of the original wire and the image wires. The results would be as shown 
in fig.1.5. None of the equipotential surfaces is a plane. The solution is to use an infinite set of equally 
spaced wires charged alternately positive and negative, as shown in Fig. 1.6. The symmetry of this set 
of wires is such that there must be equipotential planes mid-way between the wires. 
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Fig. 1.5 The field pattern around a positively charged wire flanked by a pair of negatively charged 
wires. 


Fig. 1.6 The field pattern around a set of equispaced parallel wires charged alternately positive and 
negative. 


Example 1.5 


An air-spaced coaxial line has inner and outer conductors with radii a and b respectively as shown in 


fig.1.7. Show that the breakdown voltage of the line is highest when In (a/ b) =1. 


é 


Fig. 1.7: The arrangement of an air-spaced coaxial line 
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Solution 


For most practical purposes the properties of air are indistinguishable from vacuum. From the 
symmetry of the problem we note that the electric field must everywhere be radial. The field between 
the conducting cylinders is identical to that of a long, uniform, line charge q placed along the axis of 


the system. 


To find the electric field of a line charge we apply the integral form of Gauss’ equation to a Gaussian 
surface consisting of a cylinder of unit length whose radius is r and whose ends are normal to the line 
charge as shown in fig.1.8. We note that, from considerations of symmetry, the electric field must be 
acting radially outwards and depend only on the radius r. 


Fig. 1.8 A Gaussian surface for calculating the electric field strength around a line charge. 
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Let the radial component of the electric field at radius r be E,(r). On the curved surface of the cylinder 
the radial component of the electric field is constant and the flux is thus the product of the electric 
field and the area of the curved surface. 


[[ B-dA =20rx1xE, (r) (1.7) 


Ss 
The flux of the electric field through the ends of the cylinder is zero because the electric field is 
parallel to these surfaces. 


We apply Gauss’ theorem [1.5] to find the relationship between the electric field, radius (7) and the 
unknown line charge q. Since S has unit length the total charge contained within it, which is denoted 
by the right-hand side of [1.5] is just g. Thus 


QnrE, (r= (1.8) 
€ 
which can be rearranged to give 
q 1 
Ez. (r) = = (1.9) 


2m, ¥ 


Since the electric field is inversely proportional to 7, it must be greatest when the radius is least, i.e. 


when r= a. 
1 
p= (1.10) 
27, a 


The potential difference between the cylinders is found from the electric field using [1.13] 


b 
b 
1 -V, = (r)ar= sf 2dr =n 2) (1.11) 
a 0 


2m&,/ 1 27€, 
a 


The negative sign tells us that if the charge on the inner cylinder is positive then the electrostatic 
potential of the outer cylinder is negative with respect to the inner cylinder. 


The unknown charge g can be eliminated between (1.10) and (1.11) to give the potential difference in 
terms of the maximum permitted electric field and the dimensions of the line. 


pea (413 


b a max 
a 
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The condition that the potential difference should be maximum is found by differentiating the 
potential difference with respect to the ratio of the dimensions of the conductors and setting the result 
to zero. If we set R = b/a the condition can be expressed as 


ri __In(R) él, 

“(zin(®))= wo =0 (1.13) 
or 

In(R) =In(b/a) =1 (1.14) 
Example 1.6 


An air-spaced transmission line consists of two parallel cylindrical conductors each 2 mm in diameter 
with their centres 10 mm apart as shown in fig. 1.9. Calculate the maximum potential difference 
which can be applied to the conductors assuming that the electrical breakdown strength of air 


is3 MV-m". 


Fig. 1.9 A cross-sectional view of a parallel-wire transmission line. 
Solution 
Since the diameters of the wires are small compared with their separation it is reasonable to assume 


that close to the surface of each wire the field pattern is determined almost entirely by that wire. The 
equipotential surfaces close to the wires take the form of coaxial cylinders, as may be seen in Fig. 


1.10. This is equivalent to assuming that the two wires can be represented by uniform line charges + q 
along their axes. Note that this approximation is only valid if the diameters of the wires are small 
compared with the spacing between them. 
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Fig. 1.10 The field pattern around a parallel-wire transmission line 


The electric field of either wire is then given by Equation (1.9) (for 7 =1 mm) with the appropriate 
sign for g. Since the strength of the electric field of each line charge is inversely proportional to the 
distance from the charge, the greatest electric field must occur on the plane passing through the axes 
of the two conductors. Using the notation of Fig. 1.9 and Equation (1.9) the electric field on the x axis 
between the wires is found by superimposing the fields of the two wires. 


q q 


E.= (1.15) 


2me,(4d-x) 27¢,(4d+x) 


6007 'A'G Swewsds Wy] 222] @ 
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It is easy to show that this expression is a maximum on the inner surfaces of the wires (as might be 
expected from Fig. 1.10), that is, when x = + (4d - a) . The maximum permissible charge is therefore 


given by 


Yoox = 20E,E. ao) 


1.16 
0°~ max d ( ) 


The potential at points on the x axis between the wires is found from (1.15) using [1.13] 


q q q q 1 [2dtx 
V = + dx = I + G 1.17 
(x) malteras aS * 27Ey (25) cass 


where C is a constant of integration. It is convenient to choose C = 0 so that the potential is zero at the 


origin. 


The maximum permissible potential at A is obtained by substituting the maximum charge from (1.16) 


into (1.17) and setting x = (id —a) to give 


y =p ad-4) inf —") (1.18) 


A max d a 


The potential at B is -V4 so the maximum potential difference between the wires is 2Vy. Substituting 
the numbers gives the maximum voltage between the wires as 5.9 kV. 


When the wires are not thin compared with their separation the method of solution is similar but, as 
can be seen from the equipotentials in Fig. 1.10, the equivalent line charges are no longer located at 
the centres of the wires. 


Example 1.7 


A metal sphere of radius 10 mm is placed with its centre 100 mm from a flat earthed sheet of metal. 
Assuming that the breakdown strength of air is 3 MV.m’", calculate the maximum voltage which can 
be applied to the electrode without breakdown occurring. What is then the ratio of the maximum to 
the mean surface-charge density on the sphere? 
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Solution 


This problem is solved using the same procedure as the previous one. It is necessary to assume that 
the field is equivalent to that of a pair of point charges placed at the centre of the sphere and of its 
image in the plane as shown in fig.1.11. Since the diameter of the sphere is 20% of its distance from 
the plane this assumption should not be seriously in error. We note that the attraction between the 
surface charges will ensure that the charge density and the electric field are greatest at the point on 
each sphere lying closest to the other one. 


Fig. 1.11 The arrangement of the sphere and its image in the plane. 


The first step is to use Gauss’ theorem to find the electric field at a distance r from a point charge OQ. 
The problem has spherical symmetry and therefore the electric field must be constant on the surface 
of a sphere of radius r centred on the charge and directed radially outwards. The surface area of a 


sphere of radius ris 47 r’ so that from [1.5] 


dar E,(r)== (1.19) 
€ 
so that 
E(7)=7 g 5 (1.20) 
MET 


Next we use [1.13] to find the potential at a distance r from the charge. 


V(r)=-[E,(r)ar= 2 iit a 3G (1.21) 


> ar = 
4ré,r 47é,r 


where C is a constant of integration. Now 
r=d-x (1.22) 


so that the potential due to the first charge is 
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io poo, (1.23) 


4ne,(d—x) 
Similarly the potential due to the other charge is 


A ee ee ret (1.24) 


Superimposing the potentials of the two charges gives 


v (x)= ame oe (1.25) 


Ane, d—-x d+x)_ 27E, ‘(d-x)(d+x) 


At the surface of the first sphere x =d —a and 


0) d-a 
= 1.26 ee 
“) 2mé, \ a(2d -a) eae oan) 


The electric field at the surface of the sphere is found by superimposing the fields of the two charges 
using (1.20) 
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E(d-a)= g [+ | (1.28) 


a (2d—a) 


Eliminating QO between (1.27) and (1.28) gives the relationship between the breakdown field and the 
breakdown voltage 


d-a 1 1) 
Vo, HE : 1.29 
max (5 +o ( ) 


Substituting the numerical values of the quantities we find that the maximum voltage is 28.3 kV. 


From example 1.2 we know that the maximum surface charge density is 


Ce = Eq@H mg, = 8-854 10°? x310° = 26.610 °C +m” (1.30) 


ma. 


The total charge on the sphere can be computed from (1.28) 


-1 

1 1 

Q - Bisa. 47é, 2 + 2 G31) 
a (2d -a) 


so the average charge density is 


= —] 
4 1 1 ; 
Os a Ps ds eh + 2 = ry ae I + —_—*_ Cl 32) 
Ana \a  (2d-a) (2d -a) 


Cx =| 14—" | = 1,003 (1.33) 
om (2d -a) 


Example 1.8 


An electron starts with zero velocity from a cathode which is at a potential of -10 kV and then moves 
into a region of space where the potential is zero. Find its velocity. 
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Solution 


The principle of conservation of energy requires that the sum of the kinetic energy and the potential 
energy of the electron must be constant. Thus 


1 
zm taV =0 (1.34) 
where q is the charge on the electron and V is the potential relative to the cathode. The charge to mass 


ratio of an electron g/m =—1.759x10''C.kg™' and the region of zero potential has a potential relative 
q & g p p 


to the cathode +10 kV so that, rearranging (1.34) we obtain 


v= ae = 421.759 x10" x10* =59.3x10° ms” (1.35) 
m 


Note: For accelerating voltages much above 10 kV relativistic effects become important because the 
electron velocity is comparable with the velocity of light (0.2998 x 10° ms’). It is then necessary to 
use the correct relativistic expression for the kinetic energy of the electron, but the principle of the 
calculation is unchanged. 


Example 1.9 


An electron beam originating from a cathode at a potential of -10 kV has a current of 1 A and a radius 
of 10 mm. The beam passes along the axis of an earthed conducting cylinder of radius 20 mm as 
shown in fig. 1.12. Use Gauss’ theorem to find expressions for the radial electric field within the 
cylinder, and calculate the potential on the axis of the system. 


. 


Electron 
beam 


Fig. 1.12 The arrangement of an electron beam within a concentric conducting tunnel 


Note: Electron beams like this are found in the high power microwave vacuum tubes used in 
transmitters for radar, TV broadcasting and satellite communications and for powering particle 
accelerators such as the Large Hadron Collider at CERN. 
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Solution 


The velocity of the electrons is given by(1.35). The charge per unit length in an electron beam with 
current / and electron velocity v is given by 


I 
qg=-—=-16.9x10° Cm" (1.36) 
D 


The negative sign arises because the direction of the conventional current is opposite to that of the 
electron velocity. If the radius of the beam is } and it is assumed that the current density p is uniform 
within the beam then 


p=, =-53.7x10° Cm? (1.37) 
TT 


Between the electron beam and the conducting cylinder (region 2) the problem is identical to that in 
Example 1.5 and the radial electric field is given by (1.8) 


a ee (1.38) 
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Within the electron beam (region 1) Gauss’ Theorem can be applied in exactly the same way but the 
charge enclosed in unit length of a Gaussian surface of radius 7 is now 


AC eae (1.39) 


This expression replaces g in (1.38) to give the radial field in region 1 
E,(r)=—r (1.40) 


The potential in each region is found using [1.13]. In region | the result is 


q 


Ls 
27E, 


Inr +C, (1.41) 


The value of C; is chosen by requiring V; to be zero when r = a so that 


Pa ee (1.42) 
26, a 


In region 2 we have 


V,=-F [rdr=-P 1 +0, (1.43) 
” 2E, Ey 

The value of C, is chosen by setting V2 = V; when r= 6. 

P (22.2 q a 

Vo= b =r )+ In (1.44) 
° 4e, ( ) 27, (<) 

On the axis r = 0 and 

v,(0)=2 4° +—4 in] =-362 v (1.45) 

4e, 27, 


Note: This means that the electrons on the axis have a velocity slightly less than that calculated in 
(1.35) and electron velocity increases with radius. To obtain an accurate result it would be necessary 
to re-compute the electron velocities and the charge density (which now depends on 7) to obtain 
mutually consistent values. 
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Example 1.10 


Figure 1.13 shows a simplified form for the deflection plates for a low current electron beam. Given 
that the electron beam is launched from an electrode (the cathode) at a potential of -2000V and passes 
between the deflection plates as shown, estimate the angular deflection of the beam when the 
potentials of the plates are +50 V. 


Fig. 1.13 The arrangement of a pair of electrostatic deflection plates for an electron beam. 


Note: The original use of electrostatic deflection systems in cathode ray tubes for oscilloscopes is now 
obsolete but the same system can be used in machines for electron beam lithography, electron beam 
welding and scanning electron microscopes. 


Solution 


To make the problem easier we assume that the electric field is constant everywhere between the 
plates and falls abruptly to zero at the ends. Then the field between the plates is found by dividing the 
potential difference between the plates by their separation to be E, = -5000V m''. 


Because there is no x-component of E, the axial velocity of the electrons is constant and found using 
the principle of conservation of energy as in Example 1.8. 


v, =f2nV =26.5x10° ms" (1.46) 


where 7 is the charge to mass ratio of the electron. The time taken for an electron to pass along the 
length of the plates (Z) is then 


L 
t=— =1.89 ns (1.47) 
Vv 


x 


The equation of motion in the y direction for an electron is 


=-gE (1.48) 
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where q is the magnitude of the electronic charge. The transverse acceleration of the electrons is 
constant and the y-component of velocity as they leave the plates is 


E 
pa! eee me (1.49) 


y 
m 


The angle of deflection is found from the ratio of the y and x components of the velocity 


d= stn| =3.0° (1.50) 


y 


x 


Note: It is, of course, unrealistic to assume that the field between the plates has the idealized form 
chosen above. To obtain a more accurate estimate of the deflection it would be necessary to find the 
field distribution between the plates by solving Laplace’s equation. Equation (1.48) could then be 
integrated using a more realistic expression for E,,. 
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Example 1.11 


A simple thermionic diode consists of two plane parallel electrodes: the cathode and the anode. 
Electrons are emitted from the surface of the cathode with zero velocity and accelerated towards the 
anode which is maintained at a potential V, with respect to the cathode. If the density of electrons 
between the electrodes is great enough the space charge alters the distribution of the electric field. 
Show that, in the limit of high space-charge density, the current through the diode is proportional to 


yi * and independent of the rate at which electrons are supplied by the cathode. 


Solution 
The problem as stated is a one-dimensional problem in which the electron velocity, charge density 
and potential depend only on the position x. The motion of the electrons is governed by three 


equations: the non-relativistic velocity is found from (1.46) with the difference that Vis now a 
function of x. 


x= /2nV (x) (1.51) 
The current density is related to the charge density and the velocity by 
J = px (1.52) 


The relationship between the charge density and the potential is given by the 1-dimensional form of 
Poisson’s equation [1.24] 


d°V 
= ee (1.53) 
dx é, 
Eliminating the velocity and the charge density between these equations yields 
av oe: ee 
-= y 2 (1.54) 
dx Ey 21 
This equation can be integrated by multiplying both sides by 2 (d V/ dx) to give 
1 
avy 4st 5 
: = 5 ee ae (1.55) 
X i) a/ 7] 


where C is a constant. 
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To determine C we consider the effect of the electronic space-charge on the potential as shown in fig. 
1.14. If no electrons are present in the space between the electrodes the potential varies linearly with 
position as shown by the dashed line. When electrons are emitted from the cathode (at x = 0) they are 
drawn towards the anode gaining velocity as they go. Because the electrons are negatively charged 
they depress the electrostatic potential locally as shown by the solid curve. The limit to this process 
arises when the slope of the solid curve is zero at the origin because the electric field is zero there and 
no more electrons are drawn from the cathode. The current cannot be increased beyond this limit 
given by setting C= 0 in(1.55). 


0.8 


0.6 


0.2 


Normalised potential ( V / Va ) 


0 0.2 0.4 0.6 0.8 
Normalised position (x /d ) 


Fig. 1.14 The potential distribution in a space-charge limited diode 


Equation (1.55) can then be written 


= y* (1.56) 


which can be integrated by writing 


|v 4a= at dx (1.57) 


avin 


Performing the integration we get 


4 3 
a xe (1.58) 


Now V=0 when x = 0 and therefore the constant of integration C is zero. At the anode V = V, and x = 
d so that 
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as required. 


(1.59) 


The total current flowing in the diode is obtained by multiplying the current density by the area of the 
cathode surface. This equation, known as the Child-Langmuir Law, is of fundamental importance in 


the theory of vacuum electron devices which remain the dominant technology for generation of radio 


waves at high power levels. 


Example 1.12 


Find the potential distribution between a long thin conducting strip and a surrounding rectangular 
conducting tube, as shown in fig. 1.15, when the potential difference between them is 100 V. 


4a V= 100 


2a V=0 


2a 


Fig. 1.15 The arrangement of conductors for this problem. 
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Solution 


The problem may be simplified by observing that the solution is the same in each quadrant, subject to 
appropriate reflections about the planes of symmetry. One quadrant of the diagram is redrawn on an 
enlarged scale in Fig. 1.17 with a square mesh added to it. In this example we discuss the solution by 
hand. 


To start the solution we first write down the potentials on the electrodes and estimate them at all the 
interior mesh points. An easy way to do this is to assume that the potential varies linearly with 
position. These potentials are written along-side the mesh points as shown. Next we choose a starting 
point such as 4 and work through the mesh, generating new values of the potentials with Equation 
[1.30]. 


1 
y= 7th th +N) [1.30] 


where the definitions of the potentials are as shown in fig. 1.16 


Fig. 1.16 Basis of the finite difference calculation of potential. 


As each new value is calculated it is written down and the previous estimate crossed out. Figure 1.17 
shows the results of the first pass through the mesh working along each row from right to left. Along 
the lines PO and RS we make use of the symmetry of the field to supply the potentials at the mesh 
points outside the figure (i.e. V4 = V2 on PO and V3 = V; on RS). Check the figures for yourself and 
carry the process on for one more pass through the mesh to see how the solution develops. It is not 
necessary to retain many significant figures in the early stages of the calculation because any errors 
introduced do not stop the method from converging. If we work to two significant figures we can 
avoid the use of decimal points by choosing the electrode potentials at 0 and 100 V. The final values 
of the potentials can be scaled to any other potential difference if required. 
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100 


100 


100 


Fig. 1.17. The finite difference solution for one quadrant of the problem: The initial stages. 


The process is continued until no further changes are observed in the figures to the accuracy required. 
The final result is shown in fig. 1.18. Evidently the accuracy could be improved by using a finer 


mesh. 


100 


100 


100 


100 


100 


Fig. 1.18. The finite difference solution for one quadrant of the problem: The final solution. 


Example 1.13 


Figure 1.19 shows a square coaxial arrangement of electrodes. If the potential of the inner electrode is 
5 V above that of the outer electrode estimate the maximum and minimum values of the electric field 
in the space between the electrodes. 
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10 mm 


10 mm 


Fig. 1.19. A square coaxial arrangement of electrodes. 


Solution 


The finite difference method can be used to find the fields around two-dimensional arrangements of 
electrodes on which the potentials are specified. In this example we show how the method can be 
implemented on a spreadsheet. 


A uniform square mesh is defined such that the electrodes coincide with mesh lines. The mesh 
spacing is chosen so that it is small enough to provide a reasonably detailed approximation to the 
fields whilst not being so small that the computational time is very large. 


Cells of the spreadsheet are marked out such that one cell corresponds to each mesh point. The 
symmetry of the problem can be used to reduce the number of cells required. Thus, for the geometry 
shown above it is sufficient to find the solution for one quadrant of the problem. Special care is 
needed to ensure that the correct numbers of cells are used. Remember that the cells correspond to 
intersections between mesh lines and not to the cells enclosed by them. 
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The electrode potentials are entered into the cells corresponding to the electrodes and the formula in 
Equation [1.30] is entered into all the other cells. It is convenient to take the electrode potentials as 0 
and 100 to reduce the number of digits displayed. When symmetry has been used to reduce the size of 
the problem the formulae in the cells along symmetry boundaries make use of the fact that the 
potentials on either side of the boundary are equal. 


The formulae in the cells are then applied repeatedly (a process known as iteration) until the numbers 
in the cells cease to change. To do this the calculation options of the spreadsheet must be set to 
permit iteration. It is best to set the iteration to manual and to limit the number of iterations so that 
the evolution of the solution can be observed. The final numbers in the cells are then approximations 
to the potentials at the corresponding points in space. 


From this solution the equipotential curves can be plotted and the field components can be calculated at 
any mesh point by taking the ratio of the potential difference to the mesh step. Figure 1.20 shows the 
final result obtained in this way. An active version of this figure is available for download as an EXCEL 
file. Clicking on the Potential Map tab will show you the potential map plotted using the results of the 
calculations. The electric field lines could be sketched in at right angles to the equipotential lines. 


Fig. 1.20. Finite difference calculation of the problem. Mesh step = 0.25 mm. The red and blue areas 
contain fixed potentials. The white area contains the standard formula and the green areas use 
formulae which assume symmetry at the boundaries. 
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The method can be applied to more complicated problems including those with curved electrodes 
which do not fit the mesh and three-dimensional problems. Further information can be found in the 
literature. 
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2. Dielectric materials and capacitance 


2.1 Introduction 


This chapter provides examples of the solution of problems involving dielectric materials and the 
calculation of capacitance. The methods can also be used for air-spaced and vacuum capacitors. The 
introduction of materials also makes it possible to discuss problems in the theory of semi-conductor 
devices. 


2.2 Summary of the methods available 


Note: This information is provided here for convenience. The equation numbers in the companion 
volume Electromagnetism for Electronic Engineers are indicated by square brackets. 


Symbol Signifies Units 

€ (epsilon) Permittivity F.m' 

Ep Relative permittivity Dimensionless 
D Electric flux density =C.m~ 

GC Capacitance F 

W Stored energy J 


e Relationship between permittivity and relative permittivity 

E =f, 

e Definition of electric flux density D 

D=¢E [2.4] 


e The integral form of Gauss’ theorem (all materials) 


il D-ds=[[f pav [2.5] 


e The differential form of Gauss’ theorem (all materials) 


V:-D=p [2.6] 


e Boundary conditions 
The tangential component of E is continuous at a boundary 
The normal component of D is continuous at a boundary 
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e The definition of capacitance 


Q=CV 


e The energy stored in a capacitor 


e Energy stored in an electric field 
w=~ [[J DE av 
2 


e Finite difference method 
e Estimation of capacitance using energy methods 


Example 2.1 


[2.14] 


[2.16] 


[2.18] 


A MOS transistor is essentially a parallel-plate capacitor comprising a silicon substrate, a silicon 


dioxide insulating layer, and an aluminium gate electrode as shown in fig. 2.1. The silicon dioxide has 


relative permittivity 3.85 and dielectric strength 6.0 x 10° V.m’, the insulating layer is 0.1 jum thick, 


and the area of the gate electrode is 0.02 mm’. Estimate the capacitance between the gate and the 
g p g 


substrate and the maximum voltage which can be applied to the gate electrode. 


Fig. 2.1 Arrangement of layers in a MOS transistor. 


Solution 


Aluminium gate 
electrode 


Silicon dioxide 
insulator 


Silicon substrate 


When fringing fields are ignored the capacitance of a parallel plate capacitor can be calculated by 


using (1.5) in the form 


D=o 


(2.1) 
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If the potential difference between the electrodes is V and their separation is d then 


D=£,€)E = &,& -. (2.2) 


The total charge on either plate of the capacitor is 


Q= do = AD = “ty (2.3) 


where A is the area of one plate. The capacitance of a parallel plate capacitor is therefore, from [2.14] 


_ &,64 


c = 6.8 pF (2.4) 


when the numbers given in the question are inserted. The maximum permissible voltage difference 
between the gate and the substrate is the product of the dielectric strength (breakdown field) of the 
silicon dioxide and the thickness of the insulating layer (see [1.13]). The result is 60 V. 
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Example 2.2 


Using the results of Example 2.1 calculate the maximum charge per unit area which can be induced in 
the semiconductor material. If there are 2.0 x 10'* atoms per square metre in the first layer of the 
silicon crystal, what proportion can be ionized by applying a voltage to the gate which is one sixth of 
the breakdown voltage? 


Solution 
The maximum charge per unit area is obtained from (2.3) by setting A = 1 m’ and V= 60 V 


Come = CES =0.02 Cm? (2.5) 

0.1x10 
If the applied voltage is 10 V then the number of electrons per square metre corresponding to the 
surface charge is 


n= 2 = 2m — 2 13x10! m? (2.6) 
q 


6q 


where q is the charge on an electron. If we assume that this charge is represented by ionisation of 
atoms in the first layer of the silicon substrate then, dividing n by the number of atoms per square 
metre, we find that 1.06% of them are ionised 


Example 2.3 


A variable capacitor comprised a set of fixed plates, A, and a set of moving plates, B, as shown in Fig. 
2.2. The capacitor is used to tune the frequency of a resonant circuit which varies inversely as the 
square root of the capacitance. Assuming that the effects of fringing fields can be neglected, find the 
shape which the moving plates must have if the frequency is to be proportional to the angle @ in the 
range 20-160° and 500-1500 kHz. 


Fig. 2.2 Schematic diagram of a variable capacitor. A set of moving plates B rotates within a parallel 
set of fixed plates A. 
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Solution 


This capacitor is a special example of a parallel plate capacitor. The separation between the plates is 
fixed so we know from (2.4) that the capacitance is proportional to the area of overlap between the 
plates if fringing effects are neglected. Since the frequency must be linearly related to the angle let 


f =a+b0 (2.7) 
If fis in kHz and @ is in degrees then, substituting the extreme values given 
500=a+20b and 1500=a+160b (2.8) 


The solution of the pair of simultaneous equations (2.8) is b = a/50 so (2.7) becomes 


fe (ir) (2.9) 


Now the capacitance, and therefore the overlap of the plates, is proportional to the inverse square of 
the frequency so we may write the area of overlap as 


Ay 
A(9) =» 2.10 
(1+0/50)" — 


where Ao is a constant. If the plates are moved through a small angle dé then the change in the area of 


overlap is 

dA = -51°d0 (2.11) 

so that 

_ 4A = |(-2) suai! (2.12) 
do (1+ 6/50) 


Therefore the dependence of 7 on @ which is required is 
roc(1+0/50) ” (2.13) 


where @ is in degrees. 
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Example 2.4 


Show that the capacitance per unit length between the parallel wires shown in Fig. 2.3 is given by 
ME 


< ia((ala) 1) 


if d >> a. Calculate the capacitance per unit length if d= 20 mm anda = | 


mm. 


Fig. 2.3 Cross-sectional view of a parallel-wire transmission line 
Solution 


From Example 1.6 we know that the electric field on the x-axis is given by 


q q 
E= 1.15 
“ 2ne,(4d-x) 226, (4d+x) a 


The electrostatic potential difference between the wires is found using [1.13] 


£ = 2 
ee Be ge -£n( 4-1] (2.14) 
278, (fd-x) (4d+x) 2ne, \4d+x | a 


Then from the definition of capacitance [2.14] we obtain the result required. This result is only valid if 


d >> a because we have assumed that the equipotential surfaces are circles centred on the wires. 
Substituting the dimensions given we find that 


MEQ _ 


In(19) 


C= 9.4 pF.m'! (2.15) 
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Example 2.5 


The parallel wires of Example 2.4 are placed so that each is 10 mm from a large flat sheet of metal. 
How does this affect the capacitance? Calculate the capacitance per unit length between the wires in 
the presence of the metal sheet. Represent the result by an equivalent circuit which shows the effect of 
the presence of the sheet. 


Solution 


The effect of the metal sheet on the potential distribution around the wires is found by using the 
method of images as in Example 1.3(a). Let us number the line charges as shown in fig. 2.5. 


Fig. 2.4 Arrangement of charges and image charges for example 2.5 


What do you see? 


A skateboarder? 


We see an All Options Junior Trader leaving 
the office after a successful day. 
He has the freedom to make his own decisions. 
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The potential at a point distance r from a line charge g is given by 


q 
V=- I Cc 2.16 
oe n(r) + (2.16) 


where C is a constant of integration. Summing for all the charges we find that 


ee n( 2 (2.17) 


27 A 


where the sum of the constants of integration has been set to zero so that the zero of potential is on 
the mid-planes between the wires (shown as chain-dotted lines in fig. 2.4). Provided that d >> a we 


may write for a point on the surface of wire | on a line passing through the origin 7, =a, 


h=n= (a - a/,/2) and 7; = V2d-aata point on the surface of wire 1. Thus 


y,=—4 in PS? 4 4961 (2.18) 
27 1x 27.3 274 


The potential on the surface of wire 2 is -V; so the capacitance between wires | and 2 is therefore 


Cz | =10.6 pF.m™ (2.19) 


which is an increase of 13% compared with the previous result. This result can be represented by an 
equivalent circuit which includes the parasitic capacitance between the wires and the metal sheet as 


shown in fig. 2.5. 


oT T Cc 


Fig. 2.5 Equivalent circuit for example 2.5 


The capacitance added by the plane comprises the two parasitic capacitors in series with each other 
and in parallel with the capacitance between the wires. Thus C, = 2.4 pF. 
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Example 2.6 


A capacitor consists of two parallel conducting plates whose area is large enough compared with their 
spacing for edge effects to be negligible. The sheets are 0.1 mm apart and the space between them is 
partially filled with polythene sheet of thickness 0.09 mm, relative permittivity 2.25 and breakdown 
strength 30 MV.m’' as shown in fig. 2.6. Assuming that the breakdown strength of air is 3 MV.m", 
calculate the maximum voltage which can be applied to the capacitor. What difference does it make if 
the polythene sheet completely fills the space between the plates? 


0.01 mm a Air 


Fig. 2.6 Capacitor with mixed dielectrics 
Solution 
Let the surface charge on the plates of the capacitor be o then by [2.4] and (1.5) D=o . Since the 


normal component of D is continuous at the boundary between the dielectrics it follows that D is the 
same everywhere. The electric fields in the two dielectrics are given by [2.4] as 


ee (2.20) 
€ 
and 
D 
E oly = (2.21) 
Eye: 
The potential drops across the two materials are found by multiplying the electric field by the 
thickness. The total potential drop is therefore 
t ol 
V = Eyi,| tyig + (2.22) 
E,. 


From (2.20) and (2.21) we see that the electric field in the air is greater than that in the polythene. 
Therefore the air will break down before the polythene. The maximum permissible voltage across the 
capacitor is obtained by substituting the breakdown strength of air into (2.22) to give 150 V. 


If the polythene completely fills the space between the plates then the maximum voltage is the 
product of the breakdown strength of the polythene and its thickness: 3.0 kV. This example shows the 
importance of avoiding air gaps in high voltage insulation because they seriously reduce its 
effectiveness. 
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Example 2.7 


Figure 2.7 shows a coaxial cable with two layers of dielectric material. Find an expression for the 
electric field at any point within the dielectric. 


Fig. 2.7 A coaxial cable with two concentric layers of dielectric material between the conductors 


Solution 


As in the simpler case with a single dielectric described in section 2.2 of the textbook we assume a 
charge q per unit length on the inner conductor. Applying Gauss’ theorem gives 


p -4 [2.7] 
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From this point on we must consider the two dielectric regions separately. If the radial components of 
the electric field are E; and E> in the inner and the outer region, respectively, then 


E,=—! (2.23) 
27€, 1 

and 

E,=—! (2.24) 
27€,1 


The potential difference across the inner layer is found using [1.13] 


b 


r=-| 1 Gy =| 4 inf“) (2.25) 
27E\r 27, b 


In the same way the potential difference across the outer region is 


rs 4 }o(2) (2.26) 
27 €, c 


The overall potential difference is the sum of these two potential differences 


v -£{ Lima) Lin(e) (2.27) 


ca 
2m \ & E> 


The expressions for the electric field can then be found as eliminating g between (2.23), (2.24) and 
(2.27). 


-l 

E,= Fes{ Lin(ofa + Lin(eb) (2.28) 
Er\ & &é5 

and 

E, = He| Line) + Lin(eb) (2.29) 
E,r\ & &5 
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Although the conditions for continuity of the field components have not been involved explicitly on 
this occasion, it is easy to show that they are satisfied by this solution. The symmetry of the problem 
ensures that the electric field is entirely radial, so the tangential field is zero on each side of the 
boundary. The continuity condition for the radial component of the electric flux density is satisfied 
automatically in the application of Gauss’ theorem. 


Example 2.8 
Find the capacitance per unit length of the coaxial system shown in Fig. 2.7. 
Solution 


Using (2.27) and [2.14] 


C =24{ Lin(o/a)+in(ef)| (2.30) 
E E5 
Example 2.9 


A coaxial cable is to be made with two dielectric layers as shown in Fig. 2.7. The inner layer is made 
from a high-quality but expensive material which has dielectric strength £, and the outer layer from a 
cheaper material of dielectric strength £,. Find an expression for the outer radius, 5, of the inner layer 
if the cable is to be made as cheap as possible without reducing its maximum working voltage. 


Solution 


From (2.23) and (2.24) we know that the maximum electric field in each region occurs at the inner 
surface (i.e. the minimum radius). Thus the radius of the boundary between the two regions cannot be 
reduced beyond the point at which the two regions reach their breakdown fields simultaneously. From 
(2.23) and (2.24) this requires 


_ 5 (2.31) 
E,a 


pp 


where £, and E> are the breakdown fields in the two dielectrics. Thus the minimum value of 5 is 


ae (2.32) 


Ey & 
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Example 2.10 


Calculate the stored energy per unit length in the coaxial cable shown in Fig. 2.8 when the voltage 
between the electrodes is Vo 

a) by using the expression for the energy stored in a capacitor, and 

b) by integrating the energy stored in the electric field, 


and compare the results. 


Fig. 2.8 A coaxial cable 
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Solution 


The capacitance per unit length of a cable with a single dielectric layer can be deduced from (2.30) as 


27€ 

C= 2.33 

In(b/a) — 
Then, using Equation [2.16], the stored energy per unit length is 

2 
w=+cy2 = 78% (2.34) 
2 In(b/a) 
Alternatively, the electric field at radius r is given by eliminating g between (2.23) and (2.25) 
V. 

f=— 1 — 2.35 

rIn(b/a) ee 
The stored energy in a thin cylinder of thickness dr at radius r is then 
dW =) DE der dr=—Blo _©& (2.36) 

: (In(o/a)) ” 
Integrating with respect to r from r =a tor = b gives 
2 

wa 2e% (2.37) 

In(b/a) 
as before. 
Example 2.11 


Figure 2.9 shows the cross-section of two adjacent tracks on a printed circuit board. If the tracks run 
parallel to each other for 50 mm, estimate the capacitance between them. The relative permittivity of 
the material of the board is 6.0. 


4mm 


us 
1 mm—} fe —| fe—1 mm a 


Fig. 2.9. The arrangement of conducting tracks on a printed circuit board. 
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Solution 


The exact solution to this problem would be rather difficult, but an adequate estimate can be obtained 
by elementary methods. Consider first the approximation of the system by two parallel cylindrical 
wires as shown in fig. 2.3. We should not be too far out if we take the separation of the conductors to 
be the same in each case and the diameters of the wires to be equal to the widths of the strips on the 
printed circuit board. 


From Example 2.4 the capacitance per metre of the parallel wires in free space is 


S He 2.38 
7 In((d/a)-1) eee) 


So, for the dimensions given, C = 0.7 pF. If the wires were wholly embedded in the dielectric then this 
capacitance would be multiplied by the relative permittivity to give 4.2 pF. 


Now the capacitance between the conductors can be regarded as two capacitances in parallel, one for 
the part of the diagram above the board (C)) and the other (C>) for the part within and below it. Since 
capacitances in parallel add, C, must be just half the capacitance between the wires in free space, so 
that C; = 0.35 pF. 


The value of C> must lie somewhere between this figure and half of the capacitance between the wires 
when they are wholly embedded in the dielectric because of the finite thickness of the dielectric 
board. Thus 0.35 pF < C) < 2.1 pF. Combining these estimates for C; and C) gives for the total 


capacitance 0.7 pF < C< 2.5 pF. Although this estimate is rather crude it is probably accurate enough 
for the purpose for which it is required. 


Example 2.12 


Calculate the capacitance per unit length for the coaxial system of electrodes shown in fig. 2.10 


2a 


Fig. 2.10. Rectangular coaxial electrodes. 
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Solution 


Fig. 2.11. The electrostatic potential distribution close to the outer conductor of the arrangement of 
electrodes shown in fig. 2.10. (see Fig. 1.18). 
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Figure 2.11 shows the potentials at the mesh points close to the outer electrode taken from Fig. 1.18. 
The potential difference across AB is 23.9 V so the electric field strength is 


fie ae (2.39) 


The charge density at A is therefore 


pu, ea (2.40) 
a/4 


from (1.6). If we assume that this is the mean charge density between R and S, then the total charge on 
a strip of the outer conductor, of unit length, lying between R and S is 


q= 07 = 23.96, Cm (2.41) 


The charge on each section of the outer conductor can be calculated in this way. The results of the 
calculations are shown in Fig. 2.11. Summing these charges gives the total charge on the outer 


conductor: 629 ¢, C.m”' for a potential difference of 100 V between the electrodes. The capacitance 


per unit length is therefore 6.3¢) F.m‘'. 


When this figure is compared with the exact value (5.87 € ) it is found to be in error by about 7%. 


This is accurate enough for many purposes but greater accuracy can be obtained, if required, by using 
a finer mesh. 


Example 2.13 


Use approximations to the equipotential surfaces and field lines to obtain estimates of the capacitance 
per unit length of the electrodes shown in fig. 2.10. 


Solution 


From the finite difference solution to this problem shown in Fig. 1.18 we can plot the equipotentials, 
field lines and charge distributions as shown in Fig. 2.12. 
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Fig. 2.12. The field map and charge density distributions derived from Fig. 1.18. 


If conducting sheets were placed along the equipotentials the field pattern would remain the same. It 
follows that the whole capacitance can be thought of as being made up of a set of capacitors in series, 
with each capacitor being bounded by a pair of adjacent equipotentials. This idea can be used to 
estimate the capacitance from a set of approximate equipotential surfaces. These surfaces must be 
chosen so that the electrodes form part of the set and match the boundary conditions on the electrodes. 
Figure 2.13(a) shows one possible way of satisfying this requirement. 


Fig. 2.13. (a) Approximate equipotentials and (b) approximate field lines for the problem 


Now consider the elementary capacitor which is formed by the surfaces which are x and (x + dx) from 
the inner electrode. The perimeter of this capacitor is (4a + 8x) so that its contribution to the 
reciprocal of the capacitance per unit length between the electrodes is 


Download free ebooks at bookboon.com 


52 


Please click the advert 


Electromagnetism for Electronic Engineers - Examples 2. Dielectric materials and capicitance 


1 dx 
seater 42 


When this is integrated from x = 0 to x = a the result is 


eee (2.43) 
C8) 


so that C = 7.3 &». 


Comparison between Figures 2.12 and 2.13(a) shows that the potential gradient close to the outer 
electrode is greater in the latter case. Thus the total charge associated with the approximate 
equipotentials is greater than the true value. We therefore conclude that the actual capacitance must be 
less than the figure above. 


In the second approach we imagine the charges and their associated field lines being rearranged as 
shown in Fig. 2.13(b). Because the charges are crowded more closely together than in their 
equilibrium state we expect the potential difference between the electrodes to be increased so that the 
estimate of the capacitance will be lower than the actual value. In this case we can imagine the whole 
capacitance as being made up of a number of wedge-shaped capacitors in parallel. 


| Yow’re full of energy 
and ideas. And that’s 
just what we are looking for. 


© UBS 2010. All rights reserved. 


Looking for a career where your ideas could really make a difference? UBS's 
Graduate Programme and internships are a chance for you to experience 

for yourself what it’s like to be part of a global team that rewards your input 
and believes in succeeding together. 


Wherever you are in your academic career, make your future a part of ours 
by visiting www.ubs.com/graduates. 


www.ubs.com/graduates a UB S 


Download free ebooks at bookboon.com 


53 


Electromagnetism for Electronic Engineers - Examples 2. Dielectric materials and capicitance 


Each elementary capacitor is bounded by a pair of field lines and can be thought of as a sector of a 
coaxial capacitor. The boundaries at the ends of the element are small arcs which do not coincide 
exactly with the electrode surfaces but, provided that d0 is small, the error is not very large. The 
capacitance of the element is readily deduced from (2.33) to be 


€ 


“C= Tale /n) 


do (2.44) 


where 7; and 77 are the inner and outer radii of the element, respectively. Since the charge 
distributions on the electrodes have been assumed to be uniform it follows from the geometry of the 
problem that r,/7,=2 for every element. Equation (2.40) is then readily integrated from 6 = 0 to 0 = 


z/4 to find the capacitance of one quadrant. The whole capacitance is therefore 


TE 


In(2) 


C= =4.5¢, (2.45) 


If we take the arithmetic mean of the two values of capacitance as the best estimate, we conclude that 
C=(5.9+14)e Fm" 


This figure is about 0.5% greater than the exact value (5.87 €) and is actually much closer to it than 
the result of the finite difference calculation. It should not be expected that the mean of the upper and 
lower bounds calculated in this way will always be as accurate as this. Nevertheless, the method 
usually gives accuracy which is remarkable considering the crudity of the assumptions made. 


Example 2.14 
Calculate the capacitance per unit length of the square coaxial arrangement of electrodes shown in 
fig.2.14: 

a) by making use of the results of the finite difference solution obtained in Example 1.13, and 


b) by energy methods 


10 mm 


10 mm 


Fig. 2.14. Square coaxial electrodes 
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Solution (a) 


The capacitance between the electrodes is found in exactly the same manner as described in Example 
2.12. Details of the calculation of capacitance using the spreadsheet of the finite difference solution 
are shown in fig. 2.15. An active version of this figure is available for download as an EXCEL file. 
The potential differences between the mesh points are calculated and stored in a row above the matrix 
and a column at the right-hand side. These figures are summed (note that the figure corresponding to 
the central plane must be halved to avoid counting it twice). The sum for the top edge and the right- 
hand side is 179.0 as would be expected from the symmetry of the problem. Thus the capacitance is 


C= 8x179 
100 


&=1436, Fem” (2.46) 


The details of the calculation and the potential map can be viewed by double-clicking on the figure. 
Solution (b) 


Consider an elementary capacitor in the form of a square tube of unit length, side 2x, and thickness dx. 


a2] = (2.47) 


8E ox 


Integrating from x = 3 tox =5 


In(10/6) (2.48) 


Te 
C 8é 


By the same argument as before this must be an upper bound. i.e. C < 15.7¢9 F.m" 
Now consider an elementary capacitor in the form of a wedge of angle d@. As before 


€ 


aC 7 (10/6) 


do (2.49) 


Integrating from 0 to 27 gives the lower bound. 


27 E 
Ca a 1836 2.50 
In(10/6) ° ee 


Taking the arithmetic mean of the two bounds gives the final value of the capacitance as 
C=14.0e, Fm! +12% 


Again we note that this is very close to the finite difference solution. 
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10.7 10.1 93 84 73 62 50 37 25 1.2 SUM = 179.0 


10 11 12 13 14 15 16 17 18 19 20 


Delta V 

89.3 89.9 90.7 91.6 92.7 93.8 95.0 963 97.5 98.8 12 
78.4 79.6 81.2 83.1 85.3 876 90.1 925 95.0 97.5 2.5 
67.2 69.0 71.4 744 77.7 81.3 85.0 888 925 96.3 3.7 
55.5 57.9 61.1 65.2 69.9 749 79.9 850 90.1 95.0 5.0 
43.1 45.9 50.0 555 61.8 683 749 81.3 876 93.8 6.2 
29.8 32.7 37.5 450 534 618 69.9 77.7 85.3 92.7 7.3 
33.6 45.0 555 65.2 744 831 91.46 8.4 

22.1 37.5 50.0 61.1 71.4 81.2 90.7 9.3 

176 32.7 459 579 69.0 796 89.9 10.1 

154 298 431 555 67.2 784 89.3 10.7 

143 281 41.3 538 658 774 888 11.2 

13.6 27.0 400 526 648 767 884 11.6 


Fig. 2.15. Calculation of capacitance from the finite difference solution. 
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Example 2.15 


Figure 2.16 shows the distribution of fixed electric charges in the depletion region of a p-n junction in 
equilibrium. Assuming that the silicon has a permittivity ¢, show that the potential difference between 
the two sides of the junction is given by: 


V, = 5 Not +N,d?) 


T 


where q is the magnitude of each charge. 


Number density of charges 


Fig. 2.16. The distribution of ionized donor and acceptor states in the depletion layer of an idealized 
p-n junction. 


Solution 


The depletion region is created by the transfer of mobile charges from one side of the junction to the 
other. Thus the total number of exposed ionised charge is the same on both sides of the junction and 


Nd, = Nd, (2.51) 


Consider first a plane at position z in the n region. The charge per unit area of the junction between 


this plane and z = d, is gNp (d, — z) where gq is the charge on each ionised atom.. Applying Gauss’ 


Theorem to this charge we find that 
E,(z)=-4N, (d, -z) (2.52) 


because the electric field is zero outside the junction. The potential difference across the part of the 
depletion region in the n-type material is 


n 2 
y, = Ne | (4, 2) de= Nol 2.33) 
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The potential difference across the part of the depletion region in the p-type material is found in the 
same way and the contact potential difference is then 


_ | 2 2 
V, = 5a lod +N,d,) (2.54) 


as required. 


Example 2.16 


The formula obtained in Example 2.15 is still valid when a reverse-bias voltage is applied to the 
junction if Vo is replaced by the sum of Vo and the reverse bias V. If a small a.c. signal dV is 
superimposed on the bias voltage, the transition capacitance of the junction is given by 


c, 42 


= 2.55 
roy (2.55) 


where Q is the total charge on either side of the junction. The net charges on the two sides of the 
junction are always equal to each other. Find an expression for the capacitance Cy in terms of V, Vo 
and the constants of the material. 

Solution 

The effect of the applied reverse bias voltage is to increase the thickness of the depletion region by 


exposing more fixed ionised atoms. Since the total charge on each side of the junction must be the 
same we may substitute for d, in (2.54) from (2.51) to give 


qd N 2 
ror L{19 82 |W (ay +) (2.56) 


A 


where V is the externally applied reverse bias and x is the increase in the thickness of the depletion 
region on the n-type side of the junction. Now the total charge exposed in the n-type material is 


Q=qN,A(d, +x) (2.57) 


where 4 is the area of the junction. Differentiating (2.56) and (2.57) with respect to x gives 


# A145), +x) (2.58) 
dx € Ny 

and 

WO =i d (2.59) 
dx 
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Dividing (2.59) by (2.58) yields the expression for the transition capacitance 


(2.60) 


(2.61) 


[aNy(,.Npo\ 1 
C.=éA i4 (2.62) 
" 2e “2 [V+¥, 


We note that in (2.60) the denominator is the total thickness of the depletion layer (W) so that the 
transition capacitance may be written 


C.=— (2.63) 


which is the same as the capacitance of a parallel plate capacitor (2.4) whose thickness is equal to the 
thickness of the depletion layer. 
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3. Steady electric currents 


3. Steady electric currents 


3.1 Introduction 


This chapter provides examples of the solution of problems involving resistive materials and the 
calculation of resistance. The examples include both uniform and non-uniform electric current 


densities and show how localised heating results from concentrations of current density. 


3.2 Summary of the methods available 


Note: This information is provided here for convenience. The equation numbers in the companion 
volume Electromagnetism for Electronic Engineers are indicated in square brackets. 


Symbol 
ft (mu) 

J 

o (sigma) 
p (rho) 

R 

€ 


Signifies Units 
Mobility m.V".s" 
Current density A.m” 
Conductivity S.m! 
Resistivity Q.m 
Resistance c@) 


Electromotive force V 
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Vu 


Drift velocity of charge carriers in an electric field 


Current density in a current-carrying conductor 


J =nqv, =nguE=ocE or E=pJ 
e Ohm’s law 
V= Ga =kI 
A 
e Power dissipated in a resistor 
P=VI=PR=V"/R 
e Density of power dissipation in a current-carrying conductor 
OP=E-J dv 


Principle of conservation of charge 


IJ saa=—= fff pay 


The continuity equation (differential form of [3.13]) 


Electromotive force 


[ E-dl=é 


Cc 


[3.1] 


[3.3] 


[3.6] 


(37) 


[3.11] 


[3.13] 


[3.14] 


[3.17] 
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Example 3.1 


A current of 1 A flows in a copper wire | mm in diameter. Given that there are about 10” conduction 
electrons per cubic metre in copper, calculate their rms drift velocity. If the current is the rms value of 
a 50 Hz alternating current, how far do the electrons move along the wire? What is the power 
dissipated per metre in the wire? (p = 17.7 x 10° Q.m). 


Solution 
Using [3.3] with wire of radius r = 0.5 mm carrying current / 


I 
v, = a + =79.5 um.s | (3.1) 
nq nqzr 


where / is the rms current. The instantaneous drift velocity is 


v, (t) ae) (3.2) 


nqar’ 
so the distance travelled by an electron in half a cycle is 


= aes sin0 d@= vot 2 a r6 pm (3.3) 
nqnr O+ nqnr @ 


d 


where @=a@tand w=100zs 1. 


Using [3.6] the resistance per metre of the wire is 


R=", =0.023 Qm" (3.4) 
ur 


and from [3.7] the power dissipated per metre is 
P=PR=2B mW (3.5) 
Example 3.2 


A coaxial cable has insulation made of a slightly conducting material of resistivity p as shown in fig. 
3.1. Given that the inner and outer radii of the insulation are a and b, respectively, find an expression 
for the leakage resistance per unit length between the inner and outer conductors. 
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Fig. 3.1. Arrangement of a coaxial cable. 
Solution 


From the symmetry of the cable the flow of the leakage current must be radial. Consider a thin tube of 
radius r and thickness dr as shown. The resistance of unit length of the tube to a radial current is, from 
[3.6], the product of the resistivity and the thickness of the tube in the direction of the current flow 
divided by the area of the tube normal to the current flow 


dR=—--dr (3.6) 
2ar 


WHAT'S MISSING IN THIS EQUATION? 


__ You could be one of our future talents 


= 7 - P wo 


MAERSK INTERNATIONAL TECHNOLOGY & SCIENCE PROGRAMME 


Are you about to graduate as an engineer or geoscientist? Or have you already graduated? 
If so, there may be an exciting future for you with A.P. Moller - Maersk. 


MAERSK 
www.maersk.com/mitas 


Download free ebooks at bookboon.com 


63 


Electromagnetism for Electronic Engineers - Examples 3. Steady electric currents 


The total leakage resistance per unit length is found by adding the resistances of all such elementary 
tubes connected in series 


b 
Raf (+ dr=P in(/a) 
Qn Jr 2x (3.7) 


a 


Example 3.3 


A diffused resistor in an integrated circuit is made by diffusing a thin layer of p-type impurity into an 
n-type isolation island. Assuming that the conductivity of the p-type layer varies linearly from oo at 
the surface to zero at the interface with the n-type layer, and that the interface is d below the surface, 
find an expression for the resistance between opposite edges of a square section of p-type layer whose 
transverse dimensions are much greater than its thickness. 


Solution 


Consider an elementary sheet of the p-type layer of thickness dz lying z below the surface, as shown 
in fig. 3.2. 


Fig. 3.2 Calculation of the sheet resistance of an integrated-circuit diffused resistor. 


Because the current flow is along the length of the sheet we shall want to add up the resistances of all 
such sheets in parallel and it is best to use the conductance G = 1/R in the analysis. The conductivity 
of the sheet as a function of z is given by 


o(z)=0,(1-2z/d) (3.8) 
so the conductance of the sheet is, from the inverse of [3.6] the product of the conductivity and the 


cross-sectional area normal to the current flow divided by the length in the direction of the current 
flow 


dG = ay(1-2/4)( =) ee ised 3.9) 
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The conductance of the whole layer is thus 


d 
G=0,|(1-2/d)dz = od /2 (3.10) 


0 


so the resistance in ohms is 


R= = 


ee B.11) 
G od 


Note that the resistance is independent of the size of the square chosen. This is a general result 
applying to any sheet of resistive material. For this reason the resistances of resistive sheets are 


usually quoted in ‘Ohms per square’ (ohms/(). In this example we have, in effect, assumed unit 
potential difference across the square and calculated the sum of the currents flowing in the elementary 
sheets. 


Example 3.4 


Figure 3.3 shows the arrangement of a pinch resistor in an integrated circuit. If the p regions (shown 
shaded) are earthed and the n region is always positive, the two p-n junctions are reverse biased. A 
layer of the n region adjacent to each junction is depleted of conduction electrons and is, effectively, 
an insulator. The thickness of the depletion layer is given by 


V +Vos 
V 


Dp 


t=a (3.12) 


where JV, is a constant and V is the local potential in the channel referred to S. Given that the n channel 
has width w (into the page) and its other dimensions are as shown in Fig. 3.3 find expressions for the 
current in the channel: 
a) when the current through the resistor is small, so that the voltage Vps is much less than V¢s, 
and 
b) when the current through the resistor is not small. 


Fig. 3.3 Arrangement of a pinch resistor for an integrated circuit. 
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Solution (a) 


If Vos << Ves then, to a first approximation V = 0 and the channel thickness is constant along its 
length and equal to 2(a — t). The conductance of the channel, using the inverse of [3.6], is 


2 —t 
Ge o,w(a ) = 20,,wa 1 Vos (3.13) 
L L 


where o,, is the conductivity of the n-type material. The current in the channel is then 


2 
i= aa Wm (3.14) 


If Ves = V, then t = a, the thickness of the channel is zero and the current falls to zero. The resistor is 
then said to be pinched off: 


Solution (b) 
When Vps is not small then the thickness of the channel is not constant but the continuity equation 


requires that the current is constant everywhere. Consider a small length dx of the channel which is 
distant x from S. Its contribution to the resistance of the channel is 
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“i 
iis Bc) ee (3.15) 
20,,w(a —t) 20,,wa Vp 


The potential difference across the element is 


=] 
Fe oe do (3.16) 
20,,wa Vp 


Rearranging (3.16) and integrating along the length of the channel we have 


Vos 


L 
fp fac = | fi- ViVes | ay (3.17) 
20,,wa : Vp 


0 


Performing the integrations we obtain the required expression for the current 


_ 20,,wa 2. Vong \2 
In= Vig + 1 1+ (3.18) 
L 3\V, Vc 


When Vps << Vos it can be shown that (3.18) reduces to (3.14) by expanding the term in 


(1 +Vpy¢ / Vox) using the binomial theorem and ignoring powers of (Vps/Vcs) greater than one. 


Example 3.5 


A burglar alarm system works by detecting the change in the resistance of a thin conducting film on a 
window when the window is broken. A window is typically 1 m square, and contact with the film is 
made by conducting strips at the top and bottom of the window. Estimate the percentage change in the 
resistance which would have to be detected by the electronic circuits if a burglar cut a hole just large 
enough to pass an arm through. What difference does it make where the hole it cut? 


Solution 


When the window is broken the current flux lines must deviate to pass around the hole. If the total 
current flowing across the window is held constant then the power dissipated in the resistance of the 
window increases and therefore the resistance of the window increases. Since we wish to know the 
minimum change of resistance which must be detected it is sufficient to find a lower bound to the 
resistance of the broken window. It is shown in [3.20] that this is achieved by employing approximate 
equipotentials. 
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Let us suppose that the hole made by the burglar is 100 mm square. Now divide the whole window into 
equal squares this size, as shown in fig. 3.4, and let the resistance of each square be R. The resistance of 
one horizontal strip of such squares is that of ten squares in parallel, i.e. 0.1R. The resistance of the 
whole window before it is broken is the sum of ten such strips in series, i.e. R. Note that this result is 
expected because the resistance of a square of conducting sheet is independent of its size. 


When the window is broken let us assume that the equipotentials remain as equally spaced horizontal 
straight lines. The resistance of the window is then just the sum of the resistances between the 
equipotentials. When one of the strips has a hole in it the resistance of the strip is now that of 9 
squares in parallel, i.e. R/9 = 0.11R and the resistance of the window is9x0.1R+0.11R=1.01R 
which is a lower bound. 


Fig. 3.4 Division of the window into equal squares. 

Thus, if the electronic circuits can detect a change in the resistance of the window of at least +1%, 

they will certainly be able to detect the hole made by the burglar. At this level of approximation it 

does not matter where in the window the hole is made. A higher value of the lower bound could be 


found by using a more complex approximation to the equipotentials but it is unlikely that the 
improved estimate would justify the effort involved. 


Example 3.6 


A right-angled bend is formed in a conducting bar of square section. Estimate the length of straight 
bar which would have the same resistance as the corner. 


Solution 


Let the width and thickness of the bar be a and its resistivity be p. Then the resistance of unit length of 
a straight section of the bar is 


R=-£ (3.19) 
a 


To find an upper bound for the resistance of the bend we assume that the current flux lines are as 
shown by the broken lines in fig. 3.5. 
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Fig. 3.5 Approximate current flow lines. 


Well away from the corner we expect the current density to be uniform. Now consider an elementary 
strip of the conductor bounded by the dashed lines. The resistance of such a strip is 


AR = p(2L+2x/2) (3.20) 
adx 
The current in the strip is 
dl = a (3.21) 
a 
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The power dissipated in the strip is 


7 a 


2 2L 2 
i = dPar = A 2 aaa as (3.22) 
a 


The power dissipated in the whole bar is 
2 y 2 
ea 2 [24a = or {21+%) (3.23) 
2 4 


Now the power calculated using approximate current flux lines must always be greater than or equal 


to the actual power dissipated W =J°R where R is the resistance of the bar. Thus 


R < = [au 4 (3.24) 


To find the upper bound of the resistance of the corner we subtract the resistance of the straight 
sections of bar 


R. < o(21+%2) eEP 2 PO. = Qgagl (3.25) 
4 a 4a a 


We note that this is independent of the choice of L as might be expected from an examination of fig. 
3.5. 


To seek for a lower bound to the resistance of the corner we employ approximate equipotential 
surfaces as shown in fig. 3.6. Far away from the corner these surfaces will be at right angles to the 
sides of the bar. 


Fig. 3.6 Approximate equipotential surfaces 
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Consider an element such as the one shown which can be treated as a section of a circular current 
flow. The resistance contributed by a sub-element of thickness dr at radius r is 


gee (3.26) 
adr 
If the potentials at the ends of the bar are +V then the potential difference across the element is 
(4Vd6/z) so that the power dissipated in the whole element is 
2 "2 2 
dw" = (*) | ody = 1 4a ilycin) (3.27) 
1 prado mp 

From the geometry of the figure we observe that 
Tie (3.28) 
i, L 
for all elements. The angle 6 ranges from 0 to z/2 so the resistance of the whole bar is 

16V°a : 
w= In(1+a/L) [do “Ain(14 a/L) (3.29) 

0 
This power must be greater than or equal to the actual power dissipated in the bar given by 
We (2v)°/R so that 
2 

oS (3.30) 

w" 2aln(1+a/L) 


Subtracting the resistance of the straight sections of bar we find that the resistance of the corner is 


nS ae) 2pL op a 2L (3.31) 
° ~ 2aIn(l+a/L) a? 2in(l+a/L) a , 


This function does depend on ZL so we wish to find the value of Z which maximises it. The best way to 
do this is to plot the function using a spreadsheet. It is then found that the function has a maximum 


value R, =0.477(p/a) when a/L =4.56 and this must be a lower bound for R. 


Finally we take the geometric mean of the upper and lower bounds as our best estimate of the 
resistance of the corner so that the length of straight bar equivalent to the corner is 0.61a + 28%. 
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Example 3.7 
Find the resistance between AA and BB of the conducting strip shown in Fig. 3.7. 


A 8 
3a a 


Thickness t 
Resistivity p 


Fig. 3.7 Electric field lines and equipotential surfaces for a current flowing through a conducting strip 
with a step change in width 


Solution 


The solution can be found using the finite difference method because when a steady current flows 
through a uniform resistive medium, the electrostatic potential satisfies Laplace’s equation. Figure 3.8 
shows the implementation of the finite difference method using a spreadsheet. 
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One half of the problem has been modelled and the step size in each direction has been chosen to be 
a/5. The red and blue regions are specified as fixed potentials of 100V and OV respectively and the 
white region uses the standard finite difference formula. In the green regions it is assumed that the 
equipotential lines are normal to the boundary and that the potentials at points which are mirror 
images of each other in the boundary are equal. This procedure forces the current flow to be tangential 
to the boundary. 


An active version of this figure is available for download as an EXCEL file. Tabs showing the 
calculations and the potential map can be selected. Note that the potential map is upside down 
compared with the calculations and that the bottom right-hand corner of the map lies outside the 
conducting area and should be ignored. 


0 1 2 3 4 5 6 7 8 9 10 a 12 13 14 15 16 17 18 19 20 Delta V 

3.64 
3.64 
3.65 
3.67 
3.69 
3.71 
3.73 


3.76 


3.77 


3.78 


3.79 


SUM 37.11 
Fig. 3.8 Finite difference solution for the problem shown in Fig. 3.7. 


To find the resistance between the ends of the strip we note that each mesh square has resistance 
R= p/t. We assume that the current through each square can be calculated from the mean potential 
difference across it. The best potentials to use for this purpose are those along the left-hand edge of 
the mesh because the field is most uniform there. The current through each square is obtained using 
Ohm’s law and the currents are summed across the strip. The result for the half problem is 


(2374 (3.32) 
p 


The total current is found by doubling this figure. The resistance of the strip is then the ratio of the 
potential difference (100 V) to the current: 


AP ge? 
2x371t t 


(3.33) 
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Example 3.8 


Use energy methods to obtain an approximate value for the resistance of the conducting strip shown 
in Fig. 3.7. 


Solution 


To obtain the upper bound for R we must choose approximate current flux lines. Figure 3.9 shows one 
possibility, no current flows at all in the shaded area of the diagram and the total current is assumed to 
be uniformly distributed across any plane perpendicular to the axis. It is convenient to divide the 
problem into two parts as shown and to calculate the resistances R, and R2 separately. 


Fig. 3.9 Approximate current flow lines to find an upper bound for the resistance of the strip shown in 
Fig. 3.7. 


The total width of the strip dx is found from the geometry of the problem to be 
2 
es (3.34) 


If the total current is /, then the power dissipated the strip is the product of the resistance of the strip 
and the square of the current. 


Fe ee (3.35) 
2t(6a — x) 


Integrating from x = 0 to x = 3a gives the total power dissipated in the left-hand part of the strip 


3a 


pee): Pay (3.36) 
2t (6a —x) 2t 
0 


so that 


Download free ebooks at bookboon.com 


74 


Electromagnetism for Electronic Engineers - Examples 3. Steady electric currents 


R < = a7 In2 = 1.044 (3.37) 
1 2 
I 2t t 


The resistance of the right-hand part of the strip is found by elementary methods to be 


R, =0.52 (3.38) 


Thus the upper bound of the resistance of the strip is1.54p/t. 


Fig. 3.10 Approximate equipontential surfaces to find a lower bound for the resistance of the strip 
shown in Fig. 3.7. 


To obtain the lower bound to R we assume that the equipotential surfaces are planes perpendicular to 
the axis, as shown in Fig. 3.10. The resistances R, and R> can be calculated by elementary methods: 


R= 0.758 (3.39) 
and 
R,=0.52 (3.40) 


giving the lower bound 
R > 125% (3.41) 


We note that the result found in Example 3.7 does indeed lie between these two bounds. The 
resistance of the strip is taken to be the geometric mean of the two bounds so that the best estimate of 
R using this method is 


R=142 + 11% (3.42) 


which is very close to the result obtained with greater effort by the finite difference method. 
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4. The magnetic effects of electric currents 


4.1 Introduction 


This chapter illustrates the methods of calculating the magnetic flux density produced by given 
arrangements of current-carrying conductors. The methods are appropriate for all problems in which 
no ferromagnetic materials are present. For practical purposes all other materials can be regarded as 
free space. 


4.2 Summary of the methods available 


Note: This information is provided here for convenience. The equation numbers in the companion 
volume Electromagnetism for Electronic Engineers are indicated in square brackets. 


Symbol  Signifies Units 

Ho (mu) Primary magnetic constant H.m'! 

B Magnetic flux density T, Wb.m” 
UV Magnetic scalar potential A 
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e ©The law of force between moving charges 


F= OO, ees MO@, 


ane ea (v, A(v, Af) =F, +E, 


e Magnetic flux density due to a moving charge 


e Magnetic force on a moving charge 


F,, =Q,(v, AB) 


e Magnetic flux density due to a current element 


ap =“! (a at) 
Arr” 


e The Biot-Savart Law 


= Hel dl ar 


An r 


e The magnetic circuit law in free space 
[ B-dl=1,[[J-aa 


e Magnetic scalar potential 


1 po 
UY —U, =-—| B-dl 


Mo 


B=-, grad VU =-)VU 
div(gradV) =V*U =0 
e The integral form of Gauss’ theorem 


[[ B-aa=o0 


[4.1] 


[4.3] 


[4.4] 


[4.6] 


[4.7] 


[4.12] 


[4.13] 


[4.14] 


[4.17] 


[4.15] 
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e The differential form of Gauss’ theorem 

divB=V-B=0 [4.16] 
e Force ona current element 

F =/(dl AB) [4.18] 
Example 4.1 


Figure 4.1 shows a single circular loop of wire carrying current /. Find an expression for the magnetic 
flux density at the point P on the axis of the loop. 


Fig. 4.1. Calculation of the magnetic flux density on the axis of a circular current loop. 
Solution 


Consider the flux density at P due to the current element at R. The current element has length ad@ 


and its distance from P along the line RP is r. From the triangle ORP, 
r=a/siny (4.1) 


The flux density produced by the current element is represented by the vector dB, which must be 
normal to both RP and the current element according to the rule for vector products. The vector dB, 
therefore, lies in the plane ORP at right angles to RP, as shown in Fig. 4.1. Because RP is at right 
angles to the current element it follows from equation [4.6] that 


2 


I I 
dB = "ado =" sin? y do (4.2) 
Arr Ara 


This can be resolved into its radial and axial components 
dB,=dBsiny and dB,=dBcosy (4.3) 


Integrating around the loop we find that the total z component is 
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ie! ay (4.4) 
2a 


while B,= 0, as would be expected from symmetry. 


This method can be extended to a very wide range of problems involving currents in wires. Some 
examples of its application are given in the examples in this the chapter. The solution of problems 
involving several current-carrying conductors in free space can be made easier by using the principle 
of superposition. This is possible because the flux density is proportional to the current producing it, 
as can be seen from Equation [4.7]. 


Real coils usually have many turns of wire and one should, in theory, integrate over the cross- 
sectional area of the winding. But, if the cross-section of the winding is small compared with the 
diameter of the coil, it is possible to make an approximation. The real coil is replaced by an equivalent 
current loop having a diameter equal to the mean diameter of the coil. The error in the calculated 
values of the flux density on the axis is of the order of a few per cent when this approximation is 
made. The calculation of the off-axis flux density is trickier and involves elliptic integrals. 
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Example 4.2 


A pair of circular coils have the same axis and radius a. They are separated by a distance 2d. Find the 
value of d which makes the magnetic flux density on the axis mid-way between the coils as uniform 
as possible. 


Solution 


The field of the pair of coils is found from (4.4) using the principle of superposition. If the currents 
are in the same direction in both coils then the fields are added together. The field variation along the 
axis Close to the mid-point can be expressed as a power series in z. From the symmetry of the problem 
we note that this series can only contain even powers of z. 


The problem as stated has only one variable, namely the ratio d/a. In order to make the field as 
uniform as possible we choose the value of this ratio so that the term in z’ is zero. The field expansion 
then contains only a constant term and terms in even powers of z from the fourth power upwards. This 
choice is equivalent to requiring that the second derivative of the field of each coil with respect to z 
vanishes at the mid-point. 


Differentiating (4.4) with respect to z 


ae = Hol 4 (in y) dy (4.5) 
dz 2a dy dz 

Now 

a=ztany (4.6) 


Differentiating (4.6) with respect to y gives 


p=? ted 5 (4.7) 
d cos’ W 

Substituting for z from (4.6) in (4.7) and rearranging we obtain 

ee (4.8) 
dy sin’ y 

And, substituting for (dy/dz) in (4.5) gives 

dB, Mol 4 

—+ =-——3sin" ycos 4.9 
ae y cosy (4.9) 
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Differentiating a second time with respect to z 


d*B, 3uol d (4 dy 3Hol (4,203 2 5 sin* 
= sin” yw cos — |= 4sin” y cos’ y —sin ——. 4.10 
dz 2a° ral a v) dz 2a’ ( a id v) a ( ) 


which is zero when 
tan’y =4 (4.11) 


or 
29 (4.12) 


Therefore the field midway between the coils is as uniform as possible when the distance between 
them is equal to the radius of each coil. Coils arranged in this way are known as Helmholtz coils. 
They are sometimes used to provide uniform magnetic fields for scientific experiments. 


Example 4.3 

Find an expression for the magnetic flux density at the centre of a loop of wire in the shape of a 
regular polygon with 2N sides. Check your answer by showing that it tends to the expression for the 
flux density at the centre of a circular loop when N >. 

Solution 


The magnetic flux density at a point distant R from a length of straight wire carrying current I is 


B= Myf sina 
27R 


[4.9] 


where the wire is perpendicular to the line joining it to the point and it subtends an angle 2a there. The 
polygon is made up of 2N such pieces of wire so 


a =— (4.13) 


Substituting for a in [4.9] and multiplying by the number of sides 


ee Mplsin(z /2N) on a fel. sin(z/2N) 


(4.14) 
27R 2R al2N 
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when N->0, sin(z/2N)—>2/2N and 


Bo am (4.15) 


which is the flux density at the centre of a circular loop of radius R as given by (4.4) when y =z/2. 


Example 4.4 


Figure 4.2 shows the arrangement of a short solenoid. Find expressions for the magnetic flux density 
a) at the point P, 

b) at the centre of a solenoid whose diameter is equal to its length, and 

c) at the centre of a solenoid whose diameter is small compared with its length. 


n turns per metre fength 


Fig. 4.2. A short solenoid. 
Solution 


The solutions are obtained from (4.4) using the principle of superposition. Consider an element of the 
solenoid of length dz at a distance z from P. The current flowing in the element is n/ dz so its 


contribution to the flux density at P is 


I 
dB, = So sin y dz (4.16) 
a 


Substituting for dz using (4.8) 


I I 
dB, - Hh sine y| dy = siny dy (4.17) 
a 


sin? y 


a) The flux density at P is found by integrating (4.17) with respect to y over the length of the 
solenoid 


B 
f vf 
B, - {4 siny dy =X (cos f ~cosat) (4.18) 


a 
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b) At the centre of a solenoid of length 2a B=z/4 anda =37/4. Substituting these angles into 
(4.18) 


2°" (4.19) 


a 


c) At the centre of a solenoid whose diameter is small compared with its length az and > 0. 


Substituting these angles in (4.18) 


By = font (4.20) 
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Example 4.5 


A square coil of side a carries a current J as shown in fig. 4.3. Show that the magnetic flux density at a 
point on the axis of the coil and distance z from its plane is given by 


B= 2/2 uF a 


: ae (a’ +42? \ la? +22? 


Fig. 4.3. A square coil 
Solution 
The field at P due to the current in OS is 


pe Myf sina [4.9] 


2nr 


where r is the perpendicular distance between the wire QS and the point P. From the geometry of the 
problem 


r= (4a) +2? =1 Ja’? +4z° (4.21) 


and o is the angle RPO so that 


1 
gig a 8 (4.22) 
OP (ia) ap? ¥2a2 +42? 
So that, substituting for sin a in [4.9], 
pa bol a (4.23) 


2nr 2a? +42? 
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The z components of the fields of the four sides of the coil add together whilst the components normal 
to the z axis cancel each other out. Thus the field at P is 


B_, = 4Bsin( ZOPR) = poe (4.24) 
, 
Substituting for B and r in (4.24) gives 
2 
B= fol 2a a _ 2V2 yf a (4.25) 


dar r Ja? +427 a (a’ +42°)Va" +227 


as required. This result can be checked by setting z = 0 and comparing the result with that obtained 
from (4.14) when N=2 and R =a/2. 


Example 4.6 


Figure 4.4 shows the arrangement of a set of saddle coils. Find an expression for the magnitude and 
direction of the magnetic flux density at the point P. 


Fig. 4.4. The arrangement of a pair of saddle coils. 
Solution 


The direction of the magnetic flux due to a current in a wire is given by the right hand corkscrew rule. 
If the thumb of the right hand points in the direction of the current then the curled fingers point in the 
direction of the flux. From the symmetry of the problem the flux density at P must therefore be 
directed vertically upwards. The flux density at P is found by adding together the flux densities due to 
the individual wires. 


The flux density due to one of the straight wires is given by [4.9] in which the perpendicular distance 
between P and the wire R = a and the angle a = 2/4. The angle between B and the vertical is also 


7/4 so that the contribution to the total flux density is 
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byl sin(z/4)cos(z/4) yl 
By = = 
27a Ana 


(4.26) 


The radial component of the flux density due to an element one of the curved wires is given by (4.2) 
and (4.3) 


a 
dB. = Mo" sin? ycosy dé (4.27) 
4a 


where y = 7/4. Resolving this in the vertical direction and integrating over the length of the wire 


d. B _ Ho 


-—~ | cos@ do = ee (4.28) 


4nza 2/2 V2 82a 


Finally we sum the contributions to the flux density from all the wires to give 


I I I 
B,=4 Mo rn Myf | 3 Ly (4.29) 
4nza 82a 27a 
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Example 4.7 


Figure 4.5 shows a non-magnetic ring having inner radius R; and outer radius R on which is wound a 
uniform coil of NV turns carrying current /. Find the magnetic flux density at any point within the ring. 


N turns 


Fig 4.5. A coil wound on a ring-shaped former. 
Solution 
From considerations of symmetry, the flux density within the winding must be in the tangential 


direction and vary only with radius. Applying the magnetic circuit law [4.12] to the circular path of 
radius 7 shown by the broken line in the diagram gives 


[ B-dl=277B, = 4|[ J-dA = 4)NI (4.30) 

Hence 

Zao (4.31) 
2nr 


within the winding. Any surface spanning a path of integration lying outside the winding is crossed by 
equal numbers of conductors carrying currents in each direction so that the magnetic flux density 
outside the ring is apparently zero. 


More careful thought shows that if the winding is a single layer, then current J must encircle the ring 
once in the tangential direction. The magnetic flux density is therefore given by Equation (4.4) at 
points outside the ring and remote from it. If it is necessary to produce a ring with no flux outside it 
then a winding with an even number of layers is used. Each layer is wound in the opposite direction 
around the ring from the one beneath it so that the tangential component of the current is cancelled. 


Example 4.8 


Find an expression for the magnetic flux density at the centre of a long solenoid by considering the 
limit of a toroidal coil (the coil-former is a ring with circular cross-section) whose diameter is large 
compared with the diameter of the winding. 
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Solution 


Let the number of turns per unit length of the winding be n. Then, for a toroidal winding with a total 
of N turns and an inner radius R, 


N 
n= 
27 R, 


(4.32) 


As the radius of the torus increases while the radial thickness is kept constant the inner radius, the 
mean radius and the outer radius tend to the same value, thus r > R,. Also the azimuthal component 


of the flux density in a toroidal coil becomes the axial field component in a long solenoid. 

Substituting for NV and r in (4.31) from (4.32) 

B= pin (4.33) 
which is identical to the result obtained in (4.20) by a different method. 


Example 4.9 


Two straight wires are parallel to each other and 10 mm apart. Find the force exerted per unit length 


on each wire when the wires carry currents of +1 A. 
Solution 
Let us suppose that the current in the first wire is in the x direction and that the wires lie in the x-y 


plane. To find the flux density in the region of one wire caused by the current in the other we allow a 
to tend to 2/2 in [4.9] so that 


B= ae (4.34) 
Tr 


where r is the distance between the wires. (Note: the same result can be obtain by using the magnetic 
circuit law [4.12]. The flux density is at right angles to the wires and in the right-hand corkscrew 
sense with respect to the current in the first wire. 

The force on the second wire is given by 

F =-/(dl1AB)=-I dl, B.y [4.18] 


where y is a unit vector in the y direction. Substituting for B, we obtain 


Ml? — 42x10" x1.0° 
2ar 22 x0.01 


Fo= 


y 


=20 uN (4.35) 
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Example 4.10 


An electron moving with velocity v enters a region of space in which the magnetic flux density is B at 
right angles to the velocity of the electron, as shown in Fig. 4.6. Given that the magnetic field is 
constant for a distance L in the direction of the initial motion of the electron and then falls abruptly to 
zero, find an expression for the angular deflection of the electron motion by the magnetic field. 


Solution 


The force on a moving charged particle is given by [4.4]. This force always acts at right angles to the 
direction of motion of the particle and therefore does no work on it. The magnetic force acting on a 
particle can therefore change its direction of motion but not its speed. The simplest case is that of a 
particle moving at right angles to the direction of a uniform magnetic field, as shown in Fig. 4.6. 
Because the particle is acted on by a constant force which is perpendicular to its direction of motion, it 
must move in an arc of a circle, as shown. 


Download free ebooks at bookboon.com 


89 


Electromagnetism for Electronic Engineers - Examples 4. The magnetic effects of electric currents 


| 
| 
| 
| 


Fig. 4.6. The deflection of a stream of electrons by a transverse magnetic flux. 


The radial acceleration is provided by the magnetic force, so that 


mRo =qvB =qRoB (4.36) 

whence 

o=-1B (4.37) 
m 


This frequency is known as the cyclotron frequency, from its application in the particle accelerator of 
that name. For an electron in a magnetic flux density of 0.1 Tesla the cyclotron frequency is about 2.8 
GHz. 


If the particle has an initial component of velocity parallel to the vector B, then that component is 
unaffected by the magnetic field. The motion of the particle is then helical from the combination of 


the steady motion in the direction of B with the circular motion in the direction perpendicular to it. 


Since v= Rq@ the radius of the circular path of the particle is from (4.36) 


R- a (4.38) 

But 

Rsind=L (4.39) 

so the angular deflection is given by 

9= arcsin{ am) (4.40) 
mv 
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Example 4.11 


An electric current of density J, flows through a sample of semiconductor material in the x direction 
and a uniform magnetic field B, is imposed in the y direction as shown in fig. 4.7. The resulting 
deflection of the charge carriers produces a surface charge on the sample and an electric field E, in the 
z direction which balances the magnetic force on the charges in the steady state. Find an expression 
for the electric field strength in terms of J,, B,, the charge on each charge carrier and the density of 
charge carriers. 


J, B, 


Fig. 4.7. Current density and flux density for Example 4.11 
Solution 


In equilibrium no current flows in the z direction and the net force on a charge carrier in the z 
direction must be zero. Therefore 


gE, + qv,B, =0 (4.41) 


where gq is the charge on a charge carrier (assumed positive) and v, is the mean drift velocity of the 
carriers. Now 


J, =nqv, [3.3] 


where n is the number density of the carriers, so from (4.41) 


E, =-v,B, =-=— (4.42) 


The electric field E, is associated with a potential difference between the faces of the sample and can 
be measured by a high-impedance voltmeter. As the product J,B,, is positive, the magnitude of the 


voltage is a measure of the density of charge carriers in the sample and its sign shows their sign. This 
effect, known as the Hall Effect, is used routinely to measure the properties of samples of 
semiconductor material. It is also used as a means of measuring magnetic flux density. 
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5. The magnetic effects of iron 


5.1 Introduction 


This chapter illustrates the calculation of the fields of systems of current-carrying wires in the 
presence of iron pole-pieces using the method of images and the principle of superposition. It also 
shows the application of the magnetic circuit law to problems involving soft iron and permanent 
magnets. 


5.2 Summary of the methods available 


Note: This information is provided here for convenience. The equation numbers in the companion 
volume Electromagnetism for Electronic Engineers are indicated in square brackets. 


Symbol  Signifies Units 
(mu) Permeability Hm! 
Ly Relative permeability 

jel Magnetic field Am’ 
M Magneto-motive force A 

R Reluctance A.Wb" 
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e Definition of magnetic field 


= 
i 
| 


[5.5] 


e The general form of the Magnetic Circuit Law 
f H-dl={/ J-dA [5.7] 


e Boundary conditions 
The normal component of B is continuous across a boundary 
The tangential component of H is continuous across the boundary 


e Definition of Magneto-Motive Force (MMF) in a magnetic circuit 
M = NI 
e Definition of the reluctance of a uniform section of a magnetic circuit 


L 


R=— 
HA 


e Simplified form of the magnetic circuit law applied to iron circuits 


M=RO [5.15] 


Example 5.1 


Figure 5.1 shows a long straight wire parallel to the surface of a large flat sheet of iron and at a 
distance d from it. Given that the wire carries a current /, find expressions for the x and y components 
of the magnetic flux density at the point Q(x, y). 


Fig. 5.1 The magnetic flux lines when a current-carrying conductor P lies parallel to the surface of a 
sheet of iron. 
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Fig. 5.2 The calculation of the flux density at Q for the arrangement shown in Fig. 5.1 by replacing the 
iron sheet with an image conductor at P'’. 


Solution 


Assuming that the iron can be regarded as a perfect magnetic conductor, its effect can be represented 
by an image wire at P’as shown in Fig. 5.2. Notice carefully that, in order to make the flux lines meet 
the iron surface at right angles, the current in the image wire must be in the same direction as that in 
the original wire. The flux density at O can now be found by considering each wire in turn and 
superimposing the results. Let the angle OPQ be @. Since the field around a long straight wire is 
circularly symmetrical about the wire we can apply the magnetic circuit law to give 


I 
po (5.1) 
2ar 
The field components in the coordinate directions are then 
I 
B= _ fo cos (6.2) 
: 2ar 
and 
I 
B, =-/ sing (5.3) 
: 2ar 


If the position of QO with respect of P’ is given by 6’ andr’, then the field components due to the image 


wire are 
I 
B= Jo cos 0 (5.4) 
“ Qnrr 
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and 
Bi =- He! sin 6 
* 2Qrr' 
The total field is thus 


ae Ml ( cos@" cos@ 


, 


20 r r 
and 


i= Mol ( sing’ i sin 0 
20 F ia 
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Finally, using some simple trigonometry, we obtain 


I d —d 
B, =| 27" 74 (5.8) 
2m \ x +(y+d) x +(y-d) 
and 
I 
Pe ae ee (5.9) 
2m \ x +(y+d) x°+(y-d) 


Notice that B, is zero when y = 0 as it should be. 


Strictly speaking, the imaging plane should not quite coincide with the iron surface, because iron is 
not a perfect conductor of flux, but the difference is not usually important. Because the flux density in 
iron is limited by saturation effects it is possible to regard an iron surface as a magnetic equipotential 
only if the flux density is below the saturation level. 


Example 5.2 


A solenoid of the form shown in Fig. 5.3 may be idealized as a set of circular current loops at the 
centres of the windings. If each coil has 100 turns, mean diameter 200mm, and carries a current of 1A 
calculate the maximum and minimum values of the magnetic flux density on the axis of the solenoid. 
What is the effect of changing the coil spacing to 2a? 


Fig. 5.3 The arrangement of a solenoid having a number of coils and iron pole pieces at either end. 
Solution 


The iron pole-pieces are at symmetry planes so the coils are multiply reflected in them with the result 
that the field within the solenoid is identical to that of an infinite set of coils at the same spacing. The 
field on the axis of an individual coil is found by substituting the current, the number of turns and the 
radius of the coil in (4.4) 
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B, = 100 pi 
0.2 


sin’ y (5.10) 


From the symmetry of the problem we expect that the flux density will be a maximum in the plane of 
each coil and a minimum midway between the coils. We recall that y is the angle subtended by a coil 
at the point at which the flux density is determined and superimpose the fields of all the coils at that 
point to find the total flux density. 


a) The angle subtended at the centre of a coil by one which is distant na from it is 


na 


Y,= arctan =) (5.11) 


The total flux density at the centre of a coil is therefore 


n=l 


B, =500u, f +2) "sin? “| (5.12) 
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b) The angle subtended at a point midway between two coils by one distant (n—1/2)a from it is 


5a 
y,, = arctan | (5.13) 


and the total flux density at that point is given by 


B, = 100044 > sin’ Wn (5.14) 


n=l 


The simplest way to evaluate this is to use a spreadsheet as shown in fig. 5.4. 


pi 3.141592654 In plane of coil Midway between coils 

Uo 1.25664E-06 Nn Wy Bn Wn Bn 

Constant 0.000628319 
O 1.571 0.00063 
1 1.190 0.00101 1.373 0.00118 
2 0.896 0.00060 1.030 0.00079 
3 0.695 0.00033 0.785 0.00044 
4 0.559 0.00019 0.620 0.00025 
5 0.464 0.00011 0.507 0.00014 
6 0.395 0.00007 0.427 0.00009 
7 ~~ ~=0.343 0.00005 0.367 0.00006 
8 0.303 0.00003 0.322 0.00004 
9 0.271 0.00002 0.286 0.00003 
10 0.245 0.00002 0.257 0.00002 

Total Flux Density (Tesla) 0.003056 0.003048 


Fig. 5.4. Calculation of the flux density produced by the coils 


The flux density when the coil spacing is increased to 2a (and the distance from the pole-pieces to a) 
is calculated in exactly the same manner. The maximum and minimum flux densities are 1.564mT and 
1.554mT. An active version of this figure showing the details of the calculations is available for 
download as an EXCEL file. 


Example 5.3 


Figure 5.5 shows the arrangement of a magnetic circuit. Given that the coil has 1000 turns and carries 
a current of 10 mA, and that the dimensions of the circuit are: 

L=50 mm g=2mm 

w=10mm = 10000 
calculate the flux density in the air gap. 
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Air gap 


Fig. 5.5. A simple magnetic circuit made up of an iron core with an air gap in it. The magneto-motive 
force is supplied by a winding of N turns of wire 


Solution 


The problem can be represented approximately by the equivalent circuit shown in Fig. 5.6. 


Fig. 5.6. A network representing the magnetic circuit of Fig. 5.5 including the effects of flux leakage. 
The magnetomotive force is given by 


M=NI=10A (5.15) 


The reluctances of the iron and the air gap are, from the two terms of equation [5.17] 


a, -(44]-2 (5.16) 
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R, = — oe (5.17) 
Ly(w+ g) 72 My 


To estimate the leakage reluctance across A-A' we regard the path as a large air gap. The pole-face 
areas are (40 x 10) mm’ and the gap is 40 mm. To allow for the fringing of the flux let us suppose that 
the effective cross-sectional area is double the pole-face area. Then the leakage reluctance is given by 


1 .04 

RK, =—- sn — (5.18) 
HM, 2x0.04x0.01 fo 

Some simple network analysis then yields the flux in the air gap 

D, =0.6184, Wb (5,19) 


The very crude approximation to #, can be justified by considering the effect of ignoring leakage 


altogether. In that case ® = 0.62944, Wb — an error of 2%. An even more drastic approximation is to 
regard the iron as having negligible reluctance so that ® = 0.7244, Wb— an error of about 17%. The 


best estimate of the flux density in the air gap is obtained by putting in the numerical value of [Wo and 
dividing @ by the effective area of the air gap: B, =5.4mT. 


The exact solution of this problem would require a very large computer. Fortunately such high 
accuracy is not needed in most practical cases. The worked example shows how figures which are 
accurate to a few per cent can be obtained by some quite crude approximations. Even the most 
accurate calculation made using magnetic circuit theory is likely to be in error by a few per cent 
because of the non-linear properties of the iron. 


Example 5.4 


Figure 5.7 shows a ‘pot core’ made of Ferroxcube B4. Cores like this are used to make radio- 
frequency inductors. Given that the coil has 100 turns, estimate the maximum a.c. current which can 
be passed through it without the behaviour of the inductor becoming non-linear. The initial 
magnetization curve of Ferroxcube B4 is shown in Fig. 5.8. 


14 dia. 


N 
Ns 


Coil 
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Fig. 5.7. The arrangement of a ferrite pot core for a radio frequency inductor. The dimensions are in 
millimetres. 


H (kA m-') 


Fig. 5.8. The initial magnetization curve for the ferrite material Ferroxcube B4. 
Solution 


Since the flux circulating in the core is constant the flux density is greatest where the cross-sectional 
area is least. Two parts of the magnetic circuit must be examined: 


e Axial flux density in the central core: Area = 7 x 2.5° =19.6 mm” 


e Radial flux density at the end of the central core: Area = 7x 5x 2 =31.4mm7 


Evidently the flux density is greatest in the central core. The magnetisation curve is approximately 
linear up to a flux density of 0.15 T so the maximum flux in the core is 


® =19.6x10° x 0.15 =2.95x10° Wb (5.20) 


To calculate the maximum current in the winding using [5.7] we must compute the line integral of H 
around a closed path. It is not possible to do this exactly without knowledge of the detailed fields but 
a useful approximation can be found by breaking the magnetic circuit up into a number of sections as 
shown in fig. 5.9. The Magneto-motive force needed to drive the flux through each of these sections 
can be computed using elementary methods. This approach requires the assumption that the bunching 
of the flux as it turns the corner from one section to another can be represented satisfactorily by 
including some extra length in one or more of the sections as shown. 


Fig. 5.9. The sections of the magnetic circuit. 
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We note that, in the linear region of the curve in fig. 5.8 the relative permeability can be calculated 
using [5.5] as 


BOS 
Blip 4nx1077 x1000 


(5.21) 


In region | the flux density is the maximum flux density so H = 1000 A.m" and the magneto-motive 
force is 


M, =1000x8x107° =8 A (5.22) 


In region 2 the flux density and the magnetic field decreases inversely as the radius as the flux spreads 
out radially. Thus the flux density at radius r is 


co) 0.0025 


_ 5.23 
2(7) 22x0.002 r Of) 


600T 'A'g SweIsAs WIy] JER] @ 


iB \ 
3 It’s only an 


opportunity if 


you act on it 


IKEA.SE/STUDENT 


Download free ebooks at bookboon.com 


102 


Electromagnetism for Electronic Engineers - Examples 


5. The magnetic effects of iron 


B,(r) 1 -2,95x10* 0.0025 _ 1.563 


H,(r)= a 3 
Mey Mel) 20 x0.002 rr r 


and the magneto-motive force is 
rh 7 
M, = | H,(r)dr = 1.563n{ 7) =1.609 A 


In region 3 


_ 25°). . 
H, =H, 102-72 =122.5A.m 


so that 
M, =122.5x8x10" =0.98 A 


The total magneto-motive force is therefore 


M =84+2x1.61+0.98=12.2A 


(5.24) 


(5.25) 


(5.26) 


(5.27) 


(5.28) 


Dividing by the number of turns in the coil we find that the peak current should not exceed 120 mA if 


the behaviour of the inductor is to remain linear. 


Example 5.5 


Figure 5.10 shows a magnetic recording head made of Permalloy D. Given that the magnetic material 


has a relative permeability of 500, estimate the number of turns needed to produce a flux density of 


0.05 T in the air gap from a current of 10 mA. 
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mi 


Fig. 5.10. The magnetic core for a recording head (dimensions in millimetres). 
Solution 


Magnetic recording is not now as important as it once was for video and audio recordings but it is still 
used for archiving large amounts of computer data. The magnetic field in the air-gap of the magnetic 
core magnetises a small region of magnetic tape placed against it. 


If fringing of the magnetic field in the air gap is ignored then the flux density is constant throughout 
the circuit. The magneto-motive force is the sum of the MMF’s for the core and the air gap and is 
equal to the current in the coil multiplied by the number of turns. We take the path length in the core 
to be the perimeter of a circle having the mean diameter (11 mm). The thickness of the core is not 
needed for this calculation. Thus 


3 
Va 2), ey la oigiyt ge | — = 673 turns (5.29) 
I\ My LL 0.01 500 4z x10 
Example 5.6 


A magnetic circuit like that shown in Fig. 5.5 has the following parameters 


L=50 mm g=0.2 mm w=10mm 
N= 100 turns I=10A 


Given that the core is made of a magnetic material having the initial magnetization curve shown in 
fig. 5.10, estimate the flux density in the air gap. 
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BIT) 


Fig. 5.11. The graphical solution of Example 5.6, showing that doubling the magnetomotive force in a 
magnetic circuit does not necessarily produce a doubling of the flux density. 


Solution 


We assume that for this material the initial magnetization curve is an adequate approximation to the 


hysteresis loop. Allowing for the fringing around the air gap, 


2 
Bios -(“*£) Bie =1.04B,, 
WwW 


\ 
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so that 
B 
H, =| nr— £22. |1_ _ 500-7658, (5.31) 
My J4L 
or 
1 
Buon =——(500-H, 5.32 
iron rr iron ) ( ) 


A second relationship between Bj, and Hi,o, 1s given by the initial magnetisation curve in fig. 5.11. 
The solution to this pair of non-linear simultaneous equation is found by plotting the load line (5.32) 
plotted as AB on Fig. 5.11. The working point is 

Hiron= 85 A.m', Biron = 0.55 T 


and the flux density in the air gap is 
Bair = 0.53 T 


Bair 18 less than Bio, because of the fringing field of the gap. 

If the current in the coil is doubled so that the load line is CD, then the working point is 
Hiron= 210 Am’, Biron = 1.04 T 

and 
Bar = 1.0 T 


Doubling the current in the coil does not result in a doubling of the flux density in the air gap. This 
shows clearly the non-linear behaviour of the material. 


Example 5.7 


When a magnetic circuit is being designed it is likely that the reverse problem to that in Example 5.6 
has to be solved. The flux is specified and the ampere turns to produce it have to be calculated. 
Calculate the ampere turns needed to produce a flux of 0.14 mWb in the magnetic circuit of Example 
2.0; 


5. The magnetic effects of iron 
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Solution 


The flux density in the iron is 


3 
ea (5.33) 
iron (0 01) 


The magnetic field strength in the iron is read from Fig. 5.11: Hj;o, = 660 A.m''. Neglecting fringing, 
Bur = Bizop and 


B. 
H.. =—* =1.114x10°A.m'! (5.34) 
Mo 


Applying the magnetic circuit law [5.7], the total number of ampere turns is 


N=H.,,, 1, +H, L,,, =660x0.2+1.1x10° x 0.2x10-? =355 Ampere — turns (5.35) 


iron “iron air “air 
The number of turns used would depend upon the current available. 


Example 5.8 


A magnetic circuit consists of a soft iron yoke, which may be assumed to have infinite permeability, a 
permanent magnet 60 mm long and 2400 mm’ cross-sectional area, and an air gap 5 mm long and 
3600 mm‘ in cross-sectional area. The permanent magnet is made from Columax, which has the 
demagnetisation curve shown in Fig. 5.12. Initially a keeper made of soft magnetic material is 
inserted into the air gap and the magnet magnetised to a residual flux density of 1.35 T by means of a 
coil wound on the iron yoke. 


Determine: 


a) the flux density in the air gap when the keeper is removed if fringing around the gap can be 
neglected, and 

b) the value of the flux density in the permanent magnet if the keeper is replaced in the air gap, 
assuming that the recoil permeability of Columax is 1.810. 
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Fig. 5.12. The demagnetization curve for Columax. 
Solution 


a) When the keeper is removed the flux density in the air gap is related to that in the iron by 


air — — Bigs (5.36) 
3600 


WHAT'S MISSING IN THIS EQUATION? 


__ You could be one of our future talents 


= 7 - P wo 


MAERSK INTERNATIONAL TECHNOLOGY & SCIENCE PROGRAMME 


Are you about to graduate as an engineer or geoscientist? Or have you already graduated? 
If so, there may be an exciting future for you with A.P. Moller - Maersk. 


MAERSK 


www.maersk.com/mitas 


Download free ebooks at bookboon.com 


108 


Electromagnetism for Electronic Engineers - Examples 5. The magnetic effects of iron 


No MMF is supplied by a current in a winding, therefore 


B,, 
Lin = Hron lino + elon = 0 (5.37) 


iron “iron te EA aie air iron “iron 
0 


Substituting for B,ir in (5.36) from (5.35) and for the lengths of the iron and air paths 


B. 
Bigg mu, (5.38) 
Ly 60 2400 
so that the load line is 
Bron = —2-26210° x Hiron (5.39) 


The working point is found from the intersection of this line with the demagnetisation curve as shown 
by the dashed line in fig. 5.13. The working point is 


Biron =1.12T Hiron = -49 kA.m'! 


-60 -40 -20 
H({kAm-*) 


Fig. 5.13. The graphical solution to part (a) of example 5.8. 


b) When the keeper is replaced the magnet recoils to Hiron = 0 along a line whose slope is 1.8. Thus 
the increase in Biron iS 


AB on = 1-8 X My X 49,000 = 0.11 T (5.40) 


iron 


Thus the flux density becomes 1.23 T. 
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Example 5.9 


Figure 5.14(b) shows a cross-sectional view of a cylindrical loudspeaker magnet. Estimate the flux 
density in the air gap. 


50 dia. 


Frame 


arate 
coo AA Vp — te 


BRS SSS SSA AS SAAS 


Feroba 1 


Soft iron 


{a) (bd) 


Fig. 5.14. (a) A cross-sectiona| view of a moving-coil loudspeaker. The coil is placed in the air gap of 
the magnetic circuit so that it experiences an axial force when a current is passed through it. The 
movement of the coil is transferred to the air through the motion of the paper cone. (b) An enlarged 
view of the magnet assembly. 


Solution 
To get a rough estimate of the flux density in the air gap we neglect fringing and the reluctance of the 


soft iron parts of the circuit, and we assume that the magnet is operating on its main demagnetization 
curve. From the magnetic circuit law [5.7] we know that 


Fl os = ar Hs = Pei, (5.41) 
20 


The condition for the continuity of flux gives 


A. 

Bigor, = a Big = a 7 ol Big = : Des (5.42) 
Mes (25 —15 )x 20 

But B_,, = Mp ,;, So that 

Biss = 8 gE so = —3.8 * 10° Fs (5.43) 


Figure 5.15 shows the demagnetization curve for Feroba | with the load line plotted on it. The 


working point of the magnet isH,,, = 40 kA.m'', B= 0.15 T. The flux density in the air gap is 1.0 T. 
gp g y 


iron 
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H(kA m-') 


Fig. 5.15 The demagnetization curve for the ferrite permanent magnet material Feroba 1. The energy- 
product curve is shown dotted. 


Example 5.10 


Redesign the magnetic circuit shown in Fig. 5.14(b) so that the magnetic material is used as 
efficiently as possible. 
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Solution 


Figure 5.15 shows the demagnetization curve of Feroba | with the curve of BH against B added. The 
energy product is a maximum when B;,,, = 0.11 T and Hio, = 68 kA.m’. The flux density in the air 
gap is to be 1.0 T as before, so the cross-sectional area of the magnet must be 9 times that of the air 


gap, thus 4, =9x 60a mm”. Assuming that the inner diameter of the magnet is to remain 


iron 
unchanged, it can be shown that the outer diameter must be 55 mm, to give the correct area. Similarly 
the length of the magnet must be 


Lion = Bair Li, = a ee (5.44) 
Hof iron HK, 68000 
Example 5.11 


Figure 5.16 shows a cross-sectional view of a cylindrical loudspeaker magnet. The permanent magnet 
is made from Columax (Fig. 5.12) and operates on its demagnetization curve. Neglecting fringing 
fields, leakage flux, and the reluctance of the soft iron pole pieces, estimate the flux density in the air 
gap. Calculate the optimum magnet dimensions for this flux density. If the magnet were already in 
production would you consider this design change justified? 


Fig. 5.16. The arrangement of a Columax loudspeaker magnet. Notice how the use of a magnet 
material with a higher energy product allows the permanent magnet to be put at the centre instead of 
outside as shown in Fig. 5.13(b). 
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Solution 


The solution proceeds in exactly the same manner as that of Example 5.9. The load line is 


B, =-2.765x10° H, (5.45) 


iron iron 


The working point is Biron = 1.23 T, Hiron = -43 kA.m!. Thus the flux density in the air gap is 


iron 5° x 
B ip = Bion = 1.23-=0.56T (5.46) 
A 5xllx 


The maximum BH product for Columax is about 56 kJ.m™. This can be compared with the BH 
product at the working point: 53 kJ.m”. The saving in magnet volume by redesigning the magnet is 
about 6% and therefore the change is probably not worth making. 


Example 5.12 


Figure 5.17 shows an Alnico magnet for a magnetron, which is operated at a flux density of 0.4 T. 
How much weight could be saved by replacing it with a samarium cobalt magnet? The 
demagnetization curves of the two materials are shown in Fig. 5.18. The density of Alnico is 7300 
kg.m® and that of samarium cobalt is 8100 kg.m”. 


Fig. 5.17. The arrangement of an Alnico magnet for a magnetron. This magnet is typical of those used 
earlier generations of these radar transmitter tubes (dimensions in millimetres). 


1.0 
BIT) 


~800 -600 -400 -200 0 
H(kAm-') 


Fig. 5.18. The demagnetization curves for Alnico and samarium cobalt (SmCos). 
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Solution 


The Alnico magnet is operated at B = 0.4 T, H=- 30 kA.m'. Therefore its MMF is 
H1=30x10° x2x 40x10 ° 2 =7540 A (5.47) 


The energy product of samarium cobalt is maximum when B = 0.4 T and H = -350 kA.m". The two 
magnets must supply the same MMF and, therefore, the length of the samarium cobalt magnet is 


EL oe (5.48) 


~ 350x10° 


Both magnets have the same flux and, as it happens, the same flux density, so the cross-sectional area 
of the magnet is unchanged. The weights of the two magnets are 


W., = 2(0.05? —0.03”) x 0.05 x 7300 =1.83 kg (5.49) 


W, =0.05 x 0.02 x 0.0215 x 8100 = 0.17 kg (5.50) 


Thus the change to the more modern material leads to a substantial saving in both the size and the 
weight of the magnetron. The actual saving in weight would be less than appears at first sight because 
additional soft iron would be needed to complete the magnetic circuit with the samarium cobalt 
magnet. 
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6. Electromagnetic induction 


6.1 Introduction 
This chapter shows the solution of problems involving time-varying currents using Faraday’s Law of 


Electromagnetic Induction. The problems include examples both with and without magnetic materials. 
The results of the calculations are typically self and mutual inductances. 


6.2 Summary of the methods available 


Note: This information is provided here for convenience. The equation numbers in the companion 
volume Electromagnetism for Electronic Engineers are indicated in square brackets. 


Symbol Signifies Units 
c Electro-motive force (EMF) V 

A (Lambda) Flux linkage Wb 
L Self inductance H 

M Mutual inductance H 

W Stored magnetic energy J 

w Magnetic energy density Jam 
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e Faraday’s Law of Electromagnetic Induction 


e--#4 [6.14] 
dt 
oB 
E-di=-({°.aa 6.18 
J | ot [ ] 


e Definition of flux linkage (single turn) 


A=[B-dA [6.15] 
5 


e Definition of flux linkage (NV turns) 
A=N® [6.24] 


e Definition of self inductance 
L,=— [6.20] 


e Definition of mutual inductance 


A 
1 2, 


e Energy stored in an inductor 
Ig 
We 5 Li [6.38] 


e Energy storage in a magnetic field 


w=5B-H [6.45] 
w => [[]B-Hav [6.46] 


e Energy storage in iron 


wel eds [6.68] 


B, 
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Example 6.1 

Two coaxial circular wire loops, with radii 0.01 m and 0.5 m, lie in the same plane. Estimate their 
mutual inductance. If the smaller loop is rotated about a diameter, how does the mutual inductance 
depend upon the angle between the planes of the two loops? 


Solution 


The flux density on the axis at the centre of the larger loop is given by (4.4) with w =/2. Because 


the small loop is small compared with the large one we can assume that the magnetic flux density is 
uniform everywhere within it. Thus 


Batol (6.1) 
os 


where 7; is the radius of the larger loop. Then the flux linkage to the smaller loop is 


I 
Ay, =11,B, =m bai (6.2) 
27, 
From [6.19] 
Ay a 
M =— = 4a — = 0.395 nH (6.3) 
L, 2n 


If the smaller loop is rotated about a diameter so that its plane makes and angle a with the plane of the 
larger loop then the component of flux density normal to the small loop is 


B, =B,cosa (6.4) 


The mutual inductance is proportional to the normal component of the flux and therefore varies 
as cosa. 


Example 6.2 


Figure 6.1 shows a cross-section of adjacent tracks on a printed circuit board. Given that the tracks 
run parallel to each other for 50 mm, estimate the mutual inductance between the circuit comprising A 
and B and that comprising C and D. 
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CLUE 


Fig. 6.1. A cross-sectional view of the conductors on a printed circuit board. 


Solution 


Take the origin of coordinates to be mid-way between A and B as shown. The flux density at a 
distance r from a long straight wire carrying current / is given by (4.30) 


Baie 65) 


Assume that equal and opposite currents J flow in A and B. Then, applying the magnetic circuit law to 
A and B separately and superimposing the results, the net flux density at a point lying between C and 
D and x from the origin is from (6.5). 
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The magnetic flux is normal to the plane containing the conductors. The flux linked to the second 
circuit is found by integrating B over the area of circuit C-D. Since the circuit has only one turn the 
result is 


5d 


2 
fa ee in( 5 (6.7) 
2a |\3xa\x-5d x+5d 2a 3 
2 


where L is the length over which the two circuits are coupled by the magnetic field. Putting in the 
numbers shows that the mutual inductance between the circuits is 2.9 nH. 


Whether this mutual inductance will lead to appreciable crosstalk between the circuits of the kind 
illustrated by Fig. 6.1 depends upon such factors as their impedances and signal levels, Note that in 
solving the problem we have assumed that the source currents can be regarded as line currents and 
that the dimensions of the conductors can be neglected. These assumptions are not strictly valid, but 
the result obtained by making them has an accuracy which is quite adequate for most purposes. In 
cases like this the real question is whether the mutual inductance is likely to be big enough to cause 
trouble. 


Example 6.3 


The centres of the line and neutral conductors of a 30 A, 50 Hz, ring main cable are 5 mm apart, 
Estimate the maximum e.m.f. induced in a circuit enclosing an area of 0.1 m’ which is 2 m from the 
cable. 


Solution 


The amplitudes of the currents in the two wires are +3 0V2 A and their distances from the circuit 


under investigation are 2.005 m and 1.995 m. Thus the maximum flux density at the circuit is 
obtained using (6.6) as 


B 


= 30V2( AD) agg (6.8) 
2a 1.995 2.005 
The amplitude of the EMF induced in the circuit is obtained from [6.14] 


Je |= = AB = 20 x50%0.1x21x10" =0.67 AV (6.9) 


where A is the area of the circuit. 
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Example 6.4 

Figure 6.2 shows an iron core which may be assumed to have a constant relative permeability of 1000, 
a cross-sectional area of 100 mm”, and a mean path length around the magnetic circuit of 150 mm. 
Two coils of wire are wound on the core to form a simple transformer. Coil | has 100 turns and coil 2 


1000 turns. Calculate the self-inductances of the two windings and the mutual inductance between 
them. 


Fig.6.2. The arrangement of a simple transformer 
Solution 


Assume that a current /; flows in coil 1. Applying the magnetic circuit law [5.7] the magnetic field in 


the core is 
I 
H,= ae = 667x 1, Am (6.10) 


Then the flux in the core is 
®, = Au, uy H, =84 wWb (6.11) 


so that from [6.20] the self and mutual inductances are 


N,® 

L,=——_- =8.4mH (6.12) 
[ 

M ="221 84 mil (6.13) 


1 


If now a current J, is assumed to flow in coil 2 we get 


H, = 22% 2 _ 6670x1, Am™ (6.14) 
0.15 

®, = AuH, = 840 wWb (6.15) 

so that 
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c¢) 

i= Pa =840 mH (6.16) 
I, 

and 

M= “ues = 84 mH (6.17) 
2 


We have shown, incidentally, that Equation [6.22] is valid for this case and that M? = L,L, as would 


be expected from Equation [6.28] for the perfect coupling between the coils which has been assumed 
here. 


Example 6.5 


Figure 6.3 shows the arrangement of a ferrite pot core inductor for use at radio-frequencies. Given 
that the relative permeability of the ferrite is 130, what is the inductance of the coil? If the ferrite 
begins to saturate at a flux density of 0.15 T, what is the maximum current in the coil for linear 
operation? 
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25 |dia. 


gap 


Fig. 6.3. A cross-sectional view of an inductor made with a cylindrical pot core of magnetic material. 
Solution 

The magnetic circuit in this problem is similar to that in example 5.4 and can be treated using similar 
methods. We will assume that the reluctance attributable to the bending of the flux around the corners 


can be represented by extending the lengths of the inner and outer cylindrical sections of the core by 5 
mm. The reluctance of a section of the path in which the flux density is uniform is given by 


The reluctances of the air gap and the inner and outer cylindrical regions are therefore 


R, = 0.001 =e 50.9 (6.19) 
Mo x 0.0025 Ho 
G2 0.014 oe 55 (6.20) 
Mot X130x0.0025° Ly 
: A 
R; = — - _ 0.4 (6.21) 
[oH x 130% (0.025* = 0.015") Mo 
The reluctances of the top and bottom regions, where the flux is radial, are given by 
=| *— dr = ——In| 2 |= in(3)=23 (6.22) 
uJ), 2arh 2muh \ 270 Ly X 130 x 0.005 Mo 
The total reluctance of the circuit is therefore 
51.3 
R=R, +R, +2R, + K; = (6.23) 
Ho 
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If the current in the winding is I then, from [5.14] the flux is 
we NI 100 p61 


6.24 
R 57.3 ( ) 
The self inductance is found from [6.20] as 
N® 100° 
L=——= = 0.22 mH 6.25 
rt a" os 


To find the maximum current we note that the flux density is greatest in the region where the cross- 
sectional area is least namely the inner cylindrical region. Thus the maximum flux is 


®,. =24 B 


max max 


= 2x 0.00257 x 0.15 = 2.9 u.Wb (6.26) 
And the maximum current is, from [5.14] 


RO x _ 57.3X2.9x10° 


I max 


max ~ N 100 py 


=13A (6.27) 


Example 6.6 


A toroidal transformer core is made by winding 20 turns of Permalloy D strip 10 mm wide and 0.5 
mm thick onto a 40 mm diameter former. A toroidal primary winding of 100 turns is wound 
uniformly over the core and a secondary winding of 500 turns is wound uniformly over the primary. 
Calculate the self-inductances of the two windings and the mutual inductance between them, given 
that Permalloy D has a relative permeability of 4000. 


Solution 
From the dimensions given the finished core has an inner radius of 20 mm, an outer radius of 30 mm 
and a thickness of 10 mm. If the current in the primary winding is / then the magnetic field at radius r 


in the core can be found using the magnetic circuit law [5.7] 


H(r)= ae (6.28) 


Thus the flux in the core is 


= Af B(r)dr= 2 La 2 tn 2 (6.29) 
y 21 r, r 2a 4 


where / is the thickness of the core. The self inductance is given by [6.20] 
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2 2 
i= N,® _ UN ry | _ 400044) x100 0.01 30 32 mH (6.30) 
I 20 7 20 20 


Similarly the self inductance of the secondary winding is given by 


2 2 
eee ry )_ 400044) x 500 Oot (52 )=811 mi (6.31) 
2x i; 2x 20 


And the mutual inductance is 


vu ENN 2 _ 40004 x100x 5000.01 30 ~ 162 mH (6.32) 
20 A 2m 20 
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Example 6.7 


Figure 6.4 shows an iron ring of square cross-section. An inductor is made by winding N turns of wire 
uniformly on the ring as shown. Assuming that the permeability of the iron is constant, calculate the 
energy stored in the inductor when a current | is flowing in the winding, 

a) using equation [6.38], and 

b) using equation [6.46]. 


Fig. 6.4. The arrangement of a toroidal inductor. 
Solution 
a) Using equation [6.38]: 


When a current I is flowing in the winding the magnetic field at radius r is given by (6.28) and the 
flux circulating within the iron is, from (6.29) 


Ret 
ote yl 4) (6.33) 


y R+t 
ape a of | (6.34) 


a i R4t 
geal “ul 2] (6.35) 
2: 


b) Using equation [6.46]: 
Taking a volume element in the shape of a ring of radius 7, thickness dr and width a, the stored 
energy in the element is 


dW => (2nra) H(r) dr (6.36) 
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Substituting for H(r) from (6.28) and integrating over r from (R - ta)to (R + 4a) gives 


Bae 7 272 1 
R+1 
wanual (2) gapce “nl | (6.37) 
i TT 


An R-ta 


as before. It is important to realize that this calculation has been made possible by assuming that pt 
is a constant so that the system is linear. When the non-linearity of the iron has to be taken into 
account the inductance of the coil is not a constant and the energy stored depends upon the 


previous history of the iron. 


Example 6.8 


Figure 6.5 shows the arrangement of a transformer used to measure the current in a high-voltage 
cable. Calculate the open circuit voltage induced in the toroidal winding if the material of the ring has 
a relative permeability of 130 and the current in the central conductor is 10 A rms at | kHz. Does this 
answer change if the primary conductor is not on the axis of the transformer? 


10 mm 
diameter 


50 mm 


Fig. 6.5. A toroidal coil can be used as a transducer to measure the current in a conductor which 
passes through it. 


Solution 


When the conductor passes along the axis of the ring the magnetic field is given by the magnetic 


circuit law as 


fa (6.38) 
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If the thickness of the ring (2a) is small compared with its diameter (2R) we may assume that the 
magnetic field is constant within it and equal to the value at the mean radius R. Then the rms flux in 
the ring due to the central conductor is 


2 2 
® = wHra? = ns ays — =0.41 Wb (6.39) 


The rms EMF induced in the toroidal winding is given by 


e- Son =051V (6.40) 


When the wire is not in the centre of the ring the MMF remains the same even though the magnetic 
field is no longer uniform around the path of integration. Provided that the behaviour of the material 
of the core is linear the flux and the EMF induced will also be independent of the position of the 
central wire. This principle is used in current transformers. The flux density in the core may be 
detected, alternatively, by a Hall Effect sensor (see Example 4.11) and this method is used in clip-on 
ammeters. 
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Example 6.9 


Estimate the inductance per unit length of the coaxial line shown in Fig. 6.6. 


| ee 
_ 


Fig. 6.6. A rectangular coaxial line with a strip centre conductor. 
Solution 


Upper and lower bounds to the inductance of the line can be obtained using approximate equipotential 
and flux distributions in a manner analogous to that demonstrated for capacitance in Example 2.13. 


To obtain a lower bound to the inductance we assume approximate flux lines as shown by the dashed 
line in Fig. 6.6. Assuming that H is everywhere parallel to that line and applying the magnetic circuit 
law to this path we find that 


vi 
H (x)= ae (6.41) 
The flux contained in unit length of tube of thickness dx is 
Spite" a (6.42) 
4a+8x 
so that the self-inductance of the element dx is 
dL =—“©_ x (6.43) 


~ 4a +8x 


The inductances of the elements are effectively in series with each other so the total self-inductance is 
found by integrating (6.43) from x = 0 to x =a. The result is 


cs Fn(3) = 0.144% (6.44) 
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An upper bound to the inductance is obtained by assuming an approximate current distribution. One 
possibility, as shown in fig. 6.7, is to distribute the current uniformly over the broad faces of the outer 
conductor and over the whole of the inner conductor. The approximate magnetic equipotential 
surfaces associated with this distribution of current are also shown in fig. 6.7. Now the inductance of 
the region between each pair of equipotentials is associated with part of the total current flow. The 
total inductance is, therefore, approximately that of all such elements in parallel with each other. 


Fig. 6.7. An approximate current distribution in the rectangular coaxial line with a strip centre 
conductor. 
The reciprocal of the inductance of an element subtending an angle dé is obtained from [6.37] 


1 dé 
ian Hy In(b/a) (6.45) 


Since the outer conductor is twice the width of the inner conductor it follows that b/a = 2 for all the 
elements. Equation (6.45) is integrated over all the strips to give 


L < 2m(2) = 0224 (6.46) 
a 


as the upper bound of L. It will be noticed that the triangular regions at each side of the diagram have 


been omitted from the calculations so that (6.45) is integrated over z rather than 27. The inductances 
of these regions are large because the ratio b/a > ©. Thus they have the effect of reducing the 


estimate of L slightly and so the figure given above is certainly an upper bound on L. 


The best estimate of L is given by the average of the upper and lower bounds, so 
L=(0.18+0.04) p» (6.47) 


The exact figure is 0.170 so the approximate value is in error by about 6%. 
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Example 6.10 


Show how the LCRZ analogy can be used to deduce the results of Examples 2.13 and 6.9 from each 
other. 


Solution 


In Example 2.13 it was shown that 


ee es 6.48 
m(2) in) os 
so the geometrical parameter I’ must lie in the range 

1 8 
< T < —, 6.49 
in(2) in(3) ” 
Now, using the LCRZ analogy 
| eo [6.56] 


Tr 
so the inductance must lie in the range 
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Mo Mo 
; In(3) < L < ~ in(2) (6.50) 


which agrees exactly with the results of Example 6.9. 
Example 6.11 


Estimate the inductance per unit length of the square coaxial arrangement of conductors shown in fig. 
6.8. 


Dimensians 
in millimetres 


Fig.6.8. A square coaxial arrangement of electrodes 
Solution 
The solution to this problem proceeds in exactly the same manner as that of example 6.9. 


An upper bound is obtained from (6.45) by noting that each side of the line subtends an angle of 90° 
on the axis and that the ratio of the width of the outer sides to the inner sides is 5/3. Thus 


i < 5 1n(5/3) = 0.0811 (6.51) 
TT 


A lower bound is obtained in the same manner as before. The element of inductance is 


Ho 
aL = ——_~ d. 6.52 
4(a+2x) = ( ) 


where a is the length of one side of the inner conductor. Integrating from x = 0 to x =a/3 gives 
ie 4 In(5/3) = 0.0641, (6.53) 


The inductance per unit length is therefore 


L =(0.073 + 0.008) 14 (6.54) 
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7. Transmission lines 


7. Transmission lines 


7.1 Introduction 


This chapter shows the solution of problems involving pulses and waves on transmission lines. It is 


shown how the finite velocity of propagation produces distortion of pulses and changes in the 


apparent impedance as a result of the reflection of waves. 


7.2 Summary of the methods available 


Note: This information is provided here for convenience. The equation numbers in the companion 


volume Electromagnetism for Electronic Engineers are indicated in square brackets. 


Symbol Signifies Units 
k Propagation constant m' 

A (lambda) Wavelength m 

Zo Characteristic impedance Q 

I (Gamma) Voltage reflection coefficient - 

S Voltage standing wave ratio - 


e Phase velocity of a wave 


a) 
Vv, =— 


Pk 


fa 


e Phase velocity on an LC transmission line 


e Relationship between propagation constant and wavelength 


e Definition of characteristic impedance 


_V_oL 


Z 
. k 


[7.10] 


[7.11] 


[7.12] 


[7.17] 


132 


Download free ebooks at bookboon.com 


Electromagnetism for Electronic Engineers - Examples 7. Transmission lines 


e Definition of the voltage reflection coefficient 


f= oe = 0 [7.25] 


e Time-average power flow on a transmission line 


Wall [7.27] 
re 


e Definition of voltage standing wave ratio (VSWR) 


ga Vox — [7.32] 


e Apparent impedance on a line terminated by impedance Z;, 


ae Sree tan (A/) 


= 7.42 
Zo Z)+jZ,tan() a 
e Phase velocity of a TEM wave 
1 c 
ee [7.46] 
Moe Ve, 
e Wave impedance of a TEM wave 
E 
| A em [7.51] 
i» é 
Example 7.1 


A transmission line of characteristic impedance Zp terminated by a resistor R = 4 Z, as shown in 
fig.7.1 is connected to a matched source by a switch for a short time z to produce a pulse on the line. 
If the propagation time along the line is 7, where T_T, investigate the variation of voltage with time 
at each end of the line. 
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2, A B 


(a) 


Fig. 7.1. A matched source is connected to a transmission line of impedance Zp for a short time 7. 
Solution 


When the switch is closed the uncharged line appears to the source as an impedance equal to its 
characteristic impedance. The source voltage is, therefore, divided equally between the source 


impedance and the line impedance and a voltage pulse of amplitude +V, is launched on the line. After 


time 7 this pulse reaches the other end of the line where the reflection coefficient from Equation 
[7.23] 18 


19 =F 
foe ot (7.1) 
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After time 27 the pulse reflected at B returns to A. Since 7 ~—_r ‘the switch is now open and the 
reflection coefficient is 


r,= =1 Epes 


The multiple reflection of pulses can be represented by the diagram in fig. 7.2. Each time a pulse is 
reflected its amplitude is multiplied by the appropriate reflection coefficient. 


Fig. 7.2. Time of flight diagram for the circuit in fig. 7.1. 


During each reflection of the pulse the voltage appearing across the line is equal to the sum of the 
amplitudes of the incident and reflected waves. Figure 7.3 shows the voltages across the line at A and 
Basa function of time. The result of closing the switch for a short time is to produce not just one 
pulse at B but a whole string of pulses with gradually decreasing amplitudes. If the line is short so that 
T <<, then this effect shows itself as a blurring of the edges of the pulse. Evidently care has to be 
taken with the matching of lines for the digital transmission of data if errors in the information 
received are not to occur. 
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Fig. 7.3 Voltages across the line at A and B as a function of time for example 7.1 
Example 7.2 


Sketch graphs of the variation of voltage with time at A and B in the circuit of Fig. 7.1 if the source 
and load impedances are 2Z) and the switch is closed whent = 0. 


Solution 


When the switch is first closed the amplitude of the voltage pulse is+V, . The reflection coefficients at 


the source and the load are given by [7.25] 


_ 2% y-Zy _1 (73) 


07,47, 3 


The method of solution is exactly the same as in the previous example. Figure 7.4 shows the time of 
flight diagram for this case. Because the pulse persists after it has been switched on the voltages at A 
and B accumulate. Thus 


t=f% V,=24+4=0.444 
3.9 
1 1 1 
t=2T: V,=—+—+— =0.482 (7.4) 
3.9 27 
t=3T: Vea a a dad 
3.9 27 81 
and so on. 
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Fig. 7.4. Time of flight diagram for example 7.2 


Figure 7.5 shows the variation of the voltages at A and B with time. As t —> ©0 the voltages at A and B 
tend to 0.5 which is the value which would be obtained from d.c. analysis which ignores the finite 
propagation time of the pulses. 
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V4 


0.444 


~ 
os) 
~ 
an 
= 

~~ 


Fig. 7.5. Voltages across the line at A and B as a function of time for example 7.2 


Example 7.3 


A solenoid is wound uniformly at 318 turns per metre on an insulating cylindrical former 100 mm in 
diameter and 2.5 m long. Calculate its inductance per unit length. Ten tapping points are made on the 
solenoid at regular intervals, and each is connected to earth through a 0.001 wF capacitor. The 
resulting network is a cascade of symmetrical tee sections as shown in Fig. 7.6. This line is used as 
the pulse-forming network for a high-power radar transmitter in the following manner: The line is 
charged to a potential difference of 20 kV with both ends open-circuited; a matched load is then 
connected across one end and the line is discharged through it. Assuming that the line can be treated 
as a uniform transmission line, calculate the amplitude and duration of the voltage pulse supplied to 
the load. Calculate also the current flowing in the load during the pulse, and the total pulse energy. 


c 
0.001 uF 


Fig. 7.6. One section of the pulse-forming network 
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Solution 

The flux density within a long solenoid having n turns per metre is given by 

B= pnt (4.29) 
The flux linked to unit length of the solenoid is 

A= [gn I( za’) (7.5) 


where a is the radius of the former on which it is wound. Then the self inductance per unit length is 


b= A= pun" (za*)=1mH.m (7.6) 


From the information given the shunt capacitance per unit length is 


c—10x1x10% 


= 0.004 wF.m! 7.7 
25 u (7.7) 


Thus, from [7.11] and [7.17] we know that the phase velocity and characteristic impedance are 


ee ee (7.8) 


Yb 


L 
Z,= fe = 500.2 (7.9) 


When the line is discharged it behaves as a source whose impedance is equal to the characteristic 
impedance until the discharge pulse has travelled to the other end of the line and back. This is the 
reverse of the charging sequence shown in fig.7.5 of the textbook. Thus the amplitude of the voltage 
pulse applied to the load is 10 kV and its duration is 


eae (7.10) 


The current flowing into the load during the pulse is found using Ohm’s law 


_10x10° 


I =20A (7.11) 


0 


If it is assumed that the line is loss-less then the total pulse energy is equal to the energy stored in the 
capacitors when the line is fully charged 
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1 10x 0.001%10"* x(20x10°) 
1 ee ee 
2 


=2J (7.12) 


Example 7.4 
A voltage step generator whose source impedance is 50 © is connected to an oscilloscope by | m of 50 


Q, polythene-insulated, coaxial cable. If the input impedance of the oscilloscope is 1 MQ in parallel 
with 10 pF, find the waveforms appearing at the input of the oscilloscope and at the source. 


Solution 


The relative permittivity of polythene is 2.25 so, from [7.46], the phase velocity of waves on the line 


1s 


vy = =0.2x10? ms! (7.13) 


P Je. 
E, 


and the transit time for signals is 7= 5 ns. 
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Initially the source sees the cable as a matched load so the incident wave amplitude V, = +V where 


Vs; is the amplitude of the source voltage step. The amplitude of this wave remains constant because 
the cable is matched to the source and any reflected wave is absorbed by the source. Because the input 
resistance of the oscilloscope is much greater than the impedance of the source we can make the 
approximation that the effective input resistance of the oscilloscope is infinite. 


At the input to the oscilloscope the current is related to the voltage by 


(7.14) 


where t'=t+7 . The current and voltage can be expressed in terms of the amplitudes of the incident 
and reflected waves so that 


1, -1L =C“(V,4V 7.15 
ne ie a 


Then, substituting for the amplitudes of the currents in terms of the corresponding voltages and 
rearranging (7.15) 


dvi 1 
ad ZC 


Uv) (7.16) 


because V. is constant. Thus 


dvi 
les “Fe fac (7.17) 


Performing the integrations we have 


ig 


In(V, V )=5 + const (7.18) 
0 


Now the voltage across a capacitor cannot change instantaneously because, as (7.14) shows, that 
would require an infinite current to flow. Therefore V, +V_=0 when ¢'=0 so that the constant of 


integration is—In(2V, ). Substituting this expression into (7.18) and rearranging gives the equation for 
g + g p 


the amplitude of the reflected wave as a function of time 


V_=vV, [1-20 - 4) (7.19) 
0 
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we see that as t' > 00, V_—V, so that the steady-state voltage at the oscilloscope terminals is Vs 


When the system eventually reaches equilibrium the oscilloscope input voltage will be equal to Vs 
because the input resistance of the oscilloscope is much greater than the characteristic impedance of 
the cable. The time constant of the exponential is 


T=Z,C =0.5 ns (7.20) 


Thus the wave-forms at the oscilloscope input and at the source are as shown in Fig. 7.7. 


Fig. 7.7. Waveforms: (a) at the oscilloscope, and (b) at the source 


The response of the system to short pulses could be found by superimposing positive and negative 
steps at appropriate time intervals. 


Example 7.5 


A transmission line is terminated by a resistive load which is less than the characteristic impedance 

and the VSWR on the line is 2.0. Calculate the voltage reflection coefficient of the load, the ratio of 
the load impedance to the characteristic impedance of the line, and the ratio of the reflected wave to 
the incident wave, expressing your answer in decibels. 
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Solution 
By re-arranging [7.32] we find that the voltage reflection coefficient is 


—1 
i 033 (7.21) 
S+1 
Note that because the load is resistive the voltage reflection coefficient is real. Also, because the load 
Z is real and less than Zo, [ must be negative. Similarly, rearranging [7.25] 


Z_1+r_ 231 


ee ae (7.22) 
a. 12° 43-3 


Note that this ratio is equal to 1/ S if Z is real and less than Zp. If Z is real and greater than Zp the 


magnitude of the voltage reflection coefficient is the same but the ratio of impedances is equal to S. 
The ratio of the amplitude of the reflected wave to that of the incident wave is 


20log| I |=~9.5 dB (7.23) 
This is known as the Return Loss. 


Example 7.6 

A transmission line has a characteristic impedance of 50 Q, and the wavelength on the line is 150 mm 
at a frequency of | GHz. Given that the line is 20 mm long and is terminated by a 100 Q load, find the 
impedance which would be measured at the input of the line at a frequency of 1 GHz. 

Solution 


The electrical length of the line is 


20 


kl = —— x 0.02 = 0.838 radians (7.24) 
0.15 

SO 

tan(A/) ala al (7.25) 


From Equation [7.42] the input impedance of the line is 


F 


a Z, + jZ, tan(k/) 


_ 21 + jp tan(kl) - ee 


- x50 = (37.6 - 28.1) Q 7.26 
an ( i) ea 
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Example 7.7 

A coaxial cable with a characteristic impedance of 50Q is insulated with a dielectric whose relative 
permittivity is 2.7. Calculate the wavelength of signals on the line at 500 MHz, 1 GHz and 2 GHz. 
Given that a 1 metre length of this cable is terminated in a 75Q resistor, calculate the input impedance 
of the cable at each of the three frequencies given above. 


Solution 


The phase velocity of waves on this transmission line is found from [7.46] to be 


© -182x10° ms” (7.27) 


1 


The wavelengths at the three frequencies are found from [7.10] to be 365mm, 182mm and 91 mm 
respectively. The input impedances at the three frequencies are calculated in exactly the same manner 
as in the previous example to give: (33.4-1.7j) Q; (73.7+7.3j) Q; (70.1+13.4)) Q respectively. 
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Example 7.8 


A polythene-insulated coaxial cable has a characteristic impedance of 50Q and an inner conductor 
1 mm in diameter. Calculate the diameter of the outer conductor and the wavelength on the line at a 
frequency of 1 GHz. 


Solution 


The relative permittivity of polythene is 2.25. The characteristic impedance of a coaxial cable is given 
by 


Zz — nf 2 [7.47] 
2mVN é€ a 


where a and b are the diameters of the inner and outer conductors respectively and ¢ is the permittivity 
of the dielectric separating them. Hence 


in *) = 50x 24x V2.25x [= =1.25 (7.28) 
Ho 


a 
so the diameter of the outer conductor should be 
2b=2axe'* =3.5mm (7.29) 


From [7.46] the phase velocity is 


v= =0.2x10? ms”! (7.30) 


At 1 GHz the wavelength is given by [7.10] as 


A=7=02m (7.31) 
Example 7.9 


A 50Q air-spaced coaxial cable has an outer conductor with 10 mm inside diameter. Calculate the 
diameter of the inner conductor. A quarter-wave transformer is needed to match this cable to a 75Q 
load at | GHz. Assuming that the outer diameter of the transformer is the same as that of the cable, 
calculate the dimensions of the transformer. 
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Solution 


Using [7.47] we find that 


-12 
in} =S0x 20 ce = 0.834 (7.32) 
a Ar x107 


Hence the diameter of the inner conductor is 
10xe°8* = 4.34 mm (7.33) 
From [7.43] the characteristic impedance of the quarter-wave transformer must be 


Z, =V50x 75 =61.2Q (7.34) 


Thus, from (7.32), 
inf} =0.834 52 =1.021 (7.35) 
a 50 


Hence the diameter of the inner conductor is 
10xe!! =3.6 mm (7.36) 
At | GHz the free-space wavelength is 300 mm so the transformer should be 75 mm long. 


Example 7.10 

Show that an air-spaced coaxial cable is correctly matched if a sheet of resistive card having a surface 
resistance of 377 QO per square is placed across its end with a short circuit behind it. Find the 
position of the short-circuit. 

Solution 

The resistance to radial flow of electric current of a sheet having inner and outer radii a and b and 


sheet resistance Rs is calculated by considering the resistance of a ring of radial thickness dr and 
integrating from a to b. 


b 
R= | BS gees in 2] (7.37) 
a 


277 20 


a 
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Now the characteristic impedance of an air-spaced line is given by 


z,=—-. |“ (2) (7.47] 
27 \ & a 


so, if R in (7.37) is equal to Zo in [7.47] it follows that 


=. 026770 (7.38) 


Thus the resistive card forms a shunt resistance equal to the characteristic impedance of the line. If the 
line were merely continued beyond the card then the input resistance of the continuing line would be 
in parallel with the resistance of the card and the termination would not be matched. To obtain a 
matched termination the continuation of the transmission line must present an open circuit at the plane 
of the card. Therefore, from [7.43] a short circuit must be placed a quarter-wavelength behind the 


card. Note that this arrangement will only be matched exactly at the frequency at which this condition 
holds. 
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Example 7.11 


Find the phase velocity and characteristic impedance of the air-spaced two-wire transmission line 
shown in fig.7.8. 


2mm dia. 


| 20mm | 


Fig. 7.8. Air-spaced two wire transmission line 

Solution 

All uniform air-spaced transmission lines propagate TEM waves as the velocity of light. 

To find the characteristic impedance it is sufficient to determine either the inductance per unit length 
or the capacitance per unit length by analysis of the static fields. The capacitance per unit length of a 


transmission line with the dimensions shown in fig. 7.8 was found in Example 2.4 to be 9.4 pF.m"'. 
Thus, from [7.11] and [7.17] 


hi I 
Zy =,|— =—VLC = =3559 7.39 
NG 1G ae a 


Example 7.12 
Find the phase velocity and characteristic impedance for the square coaxial transmission line shown in 


fig. 7.9 if the space between the conductors is filled with a material having a relative permittivity of 
4.0. 


2mm 


4mm 


Fig. 7.9. A square coaxial transmission line 
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Solution 


The phase velocity of TEM waves on all uniform two-wire transmission lines embedded in a uniform 
dielectric is equal to the velocity of light divided by the square root of the relative permittivity of the 
dielectric (see [7.46]). Thus 


v= ; =0.15x10? m.s™ (7.40) 


To find the characteristic impedance it is sufficient to find either the capacitance or inductance per 
unit length. The capacitance per unit length of an air-spaced line having the form shown in fig.7.9 was 
found in Example 2.12 to be 6.29&  F.m’'. Note that this figure depends only on the shapes of the 
conductors and not on their absolute dimensions. The capacitance per unit length when the space 
between the conductors is filled with a dielectric is found by replacing the primary electric constant € 
by the permittivity of the dielectric. From (7.39) 


9 vé  t fm 37 =300 (7.41) 
VC 6.296,6)¢ 6.29x2\e 6.292 
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8. Maxwell’s equations and electromagnetic 
Waves 

8.1 Introduction 

The main purpose of Chapter 8 in the textbook is to summarise the laws of electromagnetism as a 


preparation for advanced study of the subject. A few examples are introduced here to illustrate some 
of the simpler applications of those laws to electromagnetic waves in free space. 


8.2 Summary of the methods available 


Note: This information is provided here for convenience. The equation numbers in the companion 
volume Electromagnetism for Electronic Engineers are indicated in square brackets. 


Symbol Signifies Units 
S The Poynting vector W.m” 


e The integral forms of Maxwell’s equations 


Il D-ds=[[f pav [8.20] 


ff B-dA =0 [8.21] 
dD 
f H-dl=[/ [3+ 2).aa [8.22] 
6B 
pE-di=—|[—-da [8.23] 


e ©The differential forms of Maxwell’s equations 


V-D=p— [8.24] 
V-B=0 [8.25] 
Ne on es [8.26] 
Ot 
VAE= — [8.27] 
Ot 
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e The constitutive relations 


D=cE [8.29] 
J=oE [8.30] 
B=uH [8.31] 


e The phase velocity of plane waves in free space 


y =! [8.37] 


J Vv EoHo 


e The wave impedance of plane waves in free space 


fo _ [#0 3770 [8.42] 
A, € 


e The Poynting vector 
S=EAH [8.45] 
e Poynting’s theorem 


The integral of (E A H) over a closed surface is equal to the instantaneous flow of 


electromagnetic power out of the volume enclosed by that surface. 


Example 8.1 


Show that plane electromagnetic waves can propagate in a uniform dielectric material and find their 
phase velocity and wave impedance. 


Solution 


The solution follows the method used in section 8.6 of the textbook. We assume that waves propagate 
in the z-direction and that the electric and magnetic fields have only x and y components respectively. 
Then, using Maxwell’s equations, we obtain as before 


OB. 
a [8.33] 
and 
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_ OH, _ 6D 


x 8.34 
OZ ot ! 


Making use of [8.29] and [8.31] with y= 44 we find that [8.35] is replaced by 


=EMy = (8.1) 


Vv, = — = = (8.2) 


Assuming propagation as exp j(@t—kz) and substituting in [8.33] we find that 
JkE, = jouH, (8.3) 


Thus the wave impedance is 


Ey _ lo My _ 1 [& 
a Ea eg [Ho 2 Po 8.4 
H, k Ho Pp € le, e ( ) 


Example 8.2 


The international guidelines for exposure of the general public to electromagnetic radiation in the 
range 400 to 2000 MHz state that it should not exceed a power density of f / 200 W.m™ where the 
frequency f is in MHz. Find the maximum electric and magnetic field strengths permissible at 2000 
MHz. 

(See: Guidelines for Limiting Exposure to Time-Varying Electric, Magnetic, and Electro-magnetic 


Fields (Up To 300 GHz), Table 7, International Commission on Non-lonizing Radiation Protection, 
1998; available at http://www. icnirp.de/documents/emfgdl.pdf ) 


Solution 


The time-average power density in an electromagnetic wave is given by 
S =4E,H, (8.5) 


where Ej and Hp are the amplitudes of the electric and magnetic fields. 
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Substituting for Hp from [8.42] 


E, =V2S x377 =V2x10377 =87 V.m* (8.6) 


Substituting for Ep in [8.42] 


E, 87 
== _=0 23 Am (8.7) 
B77 377 


Note: This information needs careful interpretation in practical cases because radio-frequency 
antennas do not radiate power uniformly in all directions and the fields at distances closer to an 
antenna than around ten wavelengths can not be approximated by plane waves 


Example 8.3 


Calculate the wave impedances of electromagnetic waves travelling in free-space, polystyrene 


(é, = 27) and alumina (e, = 8.9). Find the time-average power density in each case if the amplitude 


of the electric field is 100 V.m". 
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Solution 


The wave impedances calculated from (8.4) are 377 Q, 229 Q and 126 Q respectively. The power 
density may be calculated using (8.5) and (8.4) 


_ £ 


(8.8) 


WwW 


where Z, is the wave impedance. Substitution of the given field strength and the wave impedances 
gives 13.3, 21.8 and 39.6 W.m*” respectively. 


Example 8.4 


The reflection and transmission of plane electromagnetic waves at plane dielectric boundaries normal 
to the direction of propagation can be treated by using transmission line theory in which the electric 
and magnetic field strengths are represented by the voltage and current and the wave impedance is 
represented by the characteristic impedance. Use this method to estimate the attenuation of an 
electromagnetic wave at 2 GHz as it passes through a brick wall 120 mm thick assuming that the 
relative permittivity is 4.0 and that losses can be neglected. 


Solution 


In air the wavelength of radio waves at 2 GHz is 150 mm and the wave impedance is 377 Q. Making 
use of (8.2) and (8.4) the wavelength in brick is 75 mm and the wave impedance is 189 Q. The 
electrical thickness of the brick is 


La (8.9) 
a 75 


The problem can therefore be represented by the equivalent transmission line network shown in fig. 
8.1 in which the source is taken to be at the left of the diagram. Waves which pass through the wall 
are assumed to be absorbed in a matched load. 


Fig. 8.1. Transmission line equivalent circuit for example 8.4 
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The impedance terminating the 189 Q line at A is 377 Q. The input impedance to the line at B can be 
found from [7.42] 


_ 3+ 189 tan (10. 1) 
~ 189+ 7377 tan (10.1) 


x189 =183- /131 (8.10) 


B 


Then the voltage reflection coefficient at B is, from [7.25] 


dee Shee LO pes (8.11) 
186 — j131+377 


For unit incident power the reflected power is 
|r =0.165 (8.12) 
Therefore the amplitude of the transmitted wave in decibels relative to the incident wave is 


4=10log( 1-|P |’) =10log (0.835) =-0.79 dB (8.13) 


Example 8.5 

A long straight cylindrical wire carries a current /. Given that the wire has radius a and resistance R 
per unit length, calculate the power dissipated per unit length both directly and by integration of the 
Poynting vector. 

Solution 


(a) By direct methods. The power dissipated per unit length is/7R. 


(b) By Poynting’s theorem. Just outside the surface of the wire, the magnetic circuit law gives 


i 
A,= (8.14) 
206 
Electric field is the voltage drop per unit length 
E,=IR (8.15) 
The Poynting vector on the surface of the wire is 
PR 
s = -E,H, = (8.16) 
27a 
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The fig.8.2 shows the directions of E, and Hy. It is clear that the Poynting vector is directed into the 


surface of the wire as represented by the 


minus sign in (8.16). The power flow into unit length of the 


wire is found by integrating the Poynting vector over a cylinder just outside the wire which has radius 


a and unit length. The total power flow into the wire is 


W =-2naS,=I°R 


(8.17) 


as before. Thus the ohmic power loss in the wire can be considered the result of a steady flow of 


energy into the wire from the surrounding electromagnetic field. The circuit and field approaches are 


alternative descriptions of the same situation. 


/ 


E, we 


Fig. 8.2. Electric and magnetic fields associated with current flow in a resistive wire 
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